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ABSTRACT 
The steam generator of a DLG-9 destroyer is investigated and simu- 
lated by means of a digital computer program which was developed for this 
thesis. The Bailey Meter Company combustion control system is analyzed by 
separate control loops. The entire steam generator and control systems are 
simulated and compared with data from the DLG-9 test boiler. An attempt to 


obtain transfer functions for the boiler dependent upon steam flow is made 


and the results are analyzed. 
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L. ¿Introduction 
l1 Discussion 

The United States Navy has long been interested in the automation of 
the steam propulsion system of a surface combatant ship. This project as 
envisioned would be directed at the automation of the main propulsion plant 
to permit centrally controlled automatic operation from full power ahead to 
full power astern, and would extend to the supplementary operating procedures 
such as lighting off, warming up, transferring the electrical load, secur- 
ing the plant and so on. The decision of such a control system must en- 
sure that the system be simple and highly reliable, provide minimum response 
time in the plant to load demand changes, and provide efficient operation 
with a minimum of personnel to operate and maintain the plant. 

The Navy would like to employ a computer control installation which 
would take over in a natural progression of functions from human operators. 
The plan would proceed from unattended operation of the propulsion plant 
over the complete range of ship's power through startup and shutdown into 
the supplementary operating procedures where the cost savings and increased 
effectiveness to be had by computer automation can be studied. 

The advantage of such a computer control system would be fully realized 
on Destroyer Type ships whose operations may be divided into three situa- 
tions: 

First, the battle situation in which not only the highest possible per- 
formance of all a ship's system is required, but also the maximum ability 
to recover from any damage or derangement short of destruction of the ship. 
This, in the case of the Destroyer, includes anti-submarine warfare in which 


the ship may proceed at flank speed to the search point, then slow down 








drastically while conducting the sonar search. 

Second, the maneuvering, station-keeping and mooring situation in 
which step changes in demand, either up or down, and in rapid succession 
would be made on the propulsion plant. In this situation and in the one 
above, the fastest possible reaction time by the control system and the 
plant are desired short of dangerous overloads on any system components. 

The third situation is that of steady steaming such as when the ship 
is transiting between port and training areas, or deploying. Normal in- 
port auxiliary steaming falls into this caetgory. Under these conditions 
optimum fuel economy is a paramount consideration, and efficiency of boiler 
operation is the controlling factor. 

In addition to the automation of propulsion control, the control com- 
puter can be applied to the supplementary functions that were mentioned 
above such as lighting off and securing, which would include shifting be- 
tween the in-port and underway steaming conditions, the transferring of 
machinery loads for purposes of equalizing operating hours or in response 
to a casualty, and the operation of the electric power generating system 
which is intimately related to the management of the propulsion plant. There 
are a number of peripheral functions which lend themselves to computer con- 
trol such as, the operation and distribution of the fresh and feed water dis- 
tilling plant, the transfer and distribution of fuel oil from storage to 
service tanks and to control list, trim, and stability, and to diagnose and 
correct engineering casualties. Damage control effectiveness could be im- 
proved by processing the information influx in the computer, which could 
supply solutions to the assessed damage such as the routing of casualty 
power, counter flooding, and stability analysis. Another function of the 


computer could be the routine maintenance programming through monitoring of 








equipment hours and scheduling of maintenance to fit the ship's operating 
schedule and internal resources, 

In the area of engineering casualty control the advantage of digital 
computer control can be demonstrated. It is conceivable that the job of 
automating the plant can be accomplished by simply mechanizing the obvi- 
ous functions of a finite number of men by analog devices located at or 
near their watch stations. However, the not-obvious functions of these 
men are the difficult ones to replace. Each man monitors many plant vari- 
ables, primarily temperatures, pressures, valve positions, and the on or 
off condition of the machinery units. Each man has stored information result- 
ing from training and experience which he calls upon to make small adjust- 
ments to manipulate variables and to react correctly to excursions of impor- 
tant variables, which comprise actual or incipient engineering casualties. 
Now if the men themselves are to be on call, but not on station, the auto- 
matic control system must be able to scan all the plant variables which the 
men monitor. It must, furthermore, be programmed to initiate immediately 
casualty control measures adequate to protect the plant from damage and to 
give on-call personnel time to respond. These functions call for an informa- 
tion handling and selective response capability which no practical analog 
device could handle. This requirement on the control system stems from 
reliability problems of the machinery systems themselves and their suscepti- 
bility to derangement due to operational stresses or in battle. It is anti- 
cipated that the reliability of a digital computer control system can be 
made to far exceed anything which can be expected of the machinery itself. 

In order to undertake this project careful analysis and a mathematical 
representation of the propulsion system and particularly of the boiler, with 
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It also would be desired to implement the automatic control system on a 
destroyer-type ship which has an existing automatic combustion control device. 
This pneumatic automatic combustion control device regulates fuel, air, and 
feed water to the boiler by monitoring outlet steam pressure and feedback 
readings of the three manipulative variables. A predetermined air-to-fuel 
ratio is set into the device and some allowance is made for the longer response 
time of the forced draft system to changes in steam demand than that of the 
fuel oil system. The control device is capable of safely controlling the 
boiler over the full range of steam demands, from stop to full power and 

back to stop. 

The Naval Boiler and Turbine Laboratory has conducted extensive tests 
on a boiler from a DLG type ship, which has an automatic combustion control 
system as described above, and through a perturbation type analysis developed 
a system of transfer functions for the pneumatic control system and transfer 
functions for the boiler at two operating points Li]. 

By employing these transfer functions, the boiler and its associated 
control systems can be simulated on a computer and further study as to the 
response of the boiler to different demands can be accomplished at a minimum 
of expense. In this thesis the air, fuel, and water systems will be analyz- 
ed separately and then combined with the known transfer functions of the 
boiler at the two operating conditions, namely that of cruising condition 
and 90 percent of full power and the plant will be run as a delta model. 

Next an attempt will be made to simulate the operation of the boiler and con- 
trol systems between these two operating points so that an accurate mathe- 
matical model of the boiler and automatic combustion control system will be 
available for later implementation into the digital computer control scheme. 

The simulation will be done on a digital computer by means of a program 
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which utilizes analog computer techniques on a digital computer. This pro- 


gram was developed by the authors for use on a Control Data Corporation 1604 
Digital Computer and is an adaptation of Program JANIS [2 ]; it will be 


referred to as Program ANALOG in this thesis. 








1.2 Results and Recommendations 

The investigation of the various control loops of the Bailey Meter 
Company combustion control system showed that the fuel flow and air flow 
control systems were represented by a true mathematical model; however, 
the water flow system as represented by NBTLL1_] was oscillatory and the 
authors feel that further study of this model of the DLG-9 boiler is neces- 
sary. The combustion air flow control system was successfully converted 
from a "delta system" at two operating points to a non linear system which 
was a good mathematical representation of the entire boiler for a full range 
of operation. 

The entire steam generator plant with the combustion control system 
was simulated at the two operating points given by NBTLL11, namely, at 
the cruising condition and at the 90 per cent full power condition. The 
plant was then perturbated with a five per cent ramp change in steam flow 
both in positive and negative directions and the results were analyzed and 
found to compare favorably with expected results with the exception of the 
oscillatory water flow system. The entire plant was also tested using the 
same boiler linear transfer functions at cruising and 90 per cent condition 
and the non linear air flow control loop. A ramp change in steam flow from 
ten to 90 per cent of full power in 23 seconds was introduced into the mathe- 
matical models and the results were compared with data received from NBTL L4_] 
for a similiar test of the actual DLG-9 test boiler. The results obtained 
by using the 90 per cent transfer functions were found to agree more favor- 
ably with the test data than those results obtained by using the cruising 
conditions transfer functions. The boiler transfer functions were then made 
functions of steam flow using the cruising conditions transfer functions for 
values of steam flow below cruising conditions, and using a straight line 
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approximation between the cruising conditions and 90 per cent conditions 
transfer functions for steam flow in between those two operating points. 

The same ten to 90 per cent steam flow ramp was introduced into this simula- 
tion, but the results did not agree with the test data as well as using the 
90 per cent transfer functions. 

It is recommended that this work be continued and if possible there 
should be an attempt made to obtain transfer functions for an operating point 
of the DLG-9 boiler at or about ten per cent of full power and at another 
operating point between cruising and 90 per cent full power. With this in- 
formation a better representation of the boiler may be possible since the 
authors believe that the values of the boiler transfer functions are some 
function of steam flow and a more complete mathematical description of the 
boiler would be available with the addition of these two operating points. 
If the above recommendation cannot be accomplished, then, by the results 
of this thesis, using the 90 per cent transfer functions will result in a 
good approximation of the boiler response. 

Using these ideas, a new non-interacting type of controller can be 
designed which employs a digital control computer which will give better 
responses in steam pressure and drum water level to rapid changes in steam 
flow, than is now possible with the existing Bailey Meter Company combus- 


tion control system. 


2. Program Analog 
2.1 Discussion 

Program Analog is a digital computer program which allows problems to 
be described in an analog computer type language. Such a language, which is 
in the form of Laplace Transforms, is convenient for simulations made by 
scientists and engineers. The program actually builds up an interconnection 
of the simulated analog computer blocks (integrators, summers, multipliers, 
etc.) from a set of subroutines corresponding to each of the blocks. In 
setting up a simulation the programmer defines the interconnection structure 
between the blocks in a separate subroutine, sets the gain and parameter ad- 
justments for each block, and sets the initial conditions for the variables. 

А flow chart of the program is shown in figure 2-1. There are about 
twenty different kinds of operational blocks. They all use subscripted 
quantities for all block parameters and variables. In this way each subroutine 
may be used as many times as desired. The only constraint is that the total 
number of blocks may not exceed 98. Similar to the analog computer, special 
purpose blocks can be constructed in the same format as the others standard 
subroutines. 

The individual performing a simulation writes subroutine DIAGRAM which 
is merely a sequence of Fortran CALL statements which specify the interconnec- 
tion of the blocks. Each block used in the simulation is given a number. 
This number is used in specifying parameters for that block. Similarly each 
variable output is numbered. This number is independent of the block number, 
but for simplicity the same number is usually chosen for each. Two blocks 
are connected together by giving the same variable number to the output of 


one block and the input of the other. 


PROGRAM 


а ANALOG 





SUBROUTINE 
ÍNITCON 


SUBROUTIN E 
DIAGRAM 


(READS DATA) 


(INTERCONNECTS 
BLOCKS) 





SUBROUTINE 
INVAL 


CCALCULATES INITIAL 
VALUES) 








DELAY MAG NITD | DIVIDE 
ÍDRELAY pee LIMITER | CURVE 


Figure 2-1 Pro; 
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As an example: 





Figure 2-2(a) Example of a block sequence 


In this sequence one would be generating the formula 
y = x + x 
or in terms of the variables used in the program 
С(3) - С) ғ с(2) 
All subroutines are called in the form (N, I, J,..) where N is the block 
number, I is the output number, and J,.. the input number(s). The calling 
sequence which would be in this subroutine DIAGRAM is: 


CALL DERIVA (2,2,1) 
CALL ADDER (3,3,1,2,0) 


C(1) is the input to the DERIVA block, block number 2. The output is C(2) 
which is now the derivative of C(1). In the ADDER block, block number 3, 
the inputs are C(l) and C(2). These two variables are added, thus producing 
variable C(3). The zero in the calling sequence indicates that there is no 


third variable to be added. 
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2.2 Digital and Analog Comparison 


Some differences exist between this method of simulation and a true 


analog computer. Listed below are some advantages of each. 


2.21 Digital Advantages 


im 


Lis 


2.22 Analog 


No voltage or time scaling is required. 

Greater precision is possible, and there is no noise. 
Nonlinearities can quite easily be simulated. 

There is practically no limit to problem size. 

Delay, multiplication, limiting, etc. can be simulated 


exactly. 


Advantages 


Actual equipment can be connected in the simulation. 

For simple problems a simple computer only is needed. 
Simulation is continuous without the need for sampling; thus 
there is no sampling error involved. 

Continuous functions, such as integration, can be simulated 
exactly. 

The entire simulation proceeds simultaneously with no closed 
loop delay. 

It is possible to watch the program run and make on the spot 
adjustments in gains, initial conditions, etc. 


Problem run time is independent of problem size. 
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2.3 Analog Subroutines 

There are sixteen types of simulations available for the programmer's 
use, which are stored in the basic program as subroutines and are called 
for and interconnected in subroutine DIAGRAM. The following is a list of 
the available subroutines, their usage, the equivalent equation, and the 
data required: 
INPUT 
Usage: CALL INPUT (N, I, J) 
Equivalent Equations: C(I) = G(N) + H(N)(t) + P(N) sinQ(N)t 
Data input required: G{N), H(N), P(N), Q(N) 
REALPL 
Usage: CALL REALPL (N,I,J) 


Equivalent Laplace equation: C(I) = C(J) * G(N) 
s + H(N) 


Date input required: G(N), H(N) Note that H(N) is set equal to zero for 
Straight integration. 

DERIVA 

Usage: CALL DERIVA (N,I,J) 

Equivalent Laplace equation: C(I) = sG(N) C(J) 

Data input required: G(N) 

GAIN 

Usage: CALL GAIN (N,I,J) 

Equivalent equation: C(I) = G(N) C(J) 

Data input required:  G(N) 

ADDER 

Usage: CALL ADDER (N, I, +J, +K, +L) 

Equivalent equation: C(I) οσο) TCR) FCL) 

Data input required: None. For adding only two variables, L must be zero. 
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CMPLXPL 
Usage: CALL CMPLXPL (N,I,J) 
Equivalent Laplace equation: C(I) = G(N) C(J) 

s* + H(N) s + P(N) 
Data input required: G(N), H(N), P(N) 
DELAY1 
Usage: CALL DELAYl  (N,I,J) М = 1,2,3 

-G(N)s 

Equivalent Laplace equation: C(I) » C(J) e 
Data input required: G(N). A variable may be delayed for a maximum of 1000 
time increments. Only three Delay 1 blocks are available. 
DELAY 
Usage: CALL DELAY (N,I,J) 


0 t <G(N) C(I) = C(J) t > G(N) 


Equivalent equation: C(I) 
Data input required:  G(N) 
LIMITER 


Usage: CALL LIMITER (N,I,J) 


Equivalent equations: απ) αν «ος AO (N) 
C(I) = H(N) CON). CG). > НОМ) 
C(I) = P(N) G(N)C(J) “ғ P(N) 


Data input required: G(N), H(N), P(N) 
IDRELAY 
Usage: CALL IDRELAY (N,I,J) 


αλ ο 0. 


Equivalent equations: есі) 


C(I) H(N) СС 0 
Data input required: G(N), H(N) 
MULTPLY 


Usage: CALL MULTPLY (N,1,J,K) 
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Equivalent equation: C(I) » C(J) * C(K) 
Data input required: None. 

MAGNITD 

Usage: CALL MAGNITD (N, I, J) 

Equivalent equation: C(I) = Ic c» 

Data input required: None. 

DIVIDE 

Usage: CALL DIVIDE (N, I, J, K) 


Equivalent equation: C(J) 


= C(J 
C(K) 
Data input required: None. | 
RECIP 
Usage CALL RECIP (N, I, J) 
Equivalent equation: C(I) = G(N) 
C(J) 
Data input required:  G(N) 
ERROR 
Usage: CALL ERROR (N,I,J,K) 
Е 2 
Equivalent equation: C(I) = 1 (C(J) - C(K)) dt 
T | 
ο 


Data input required: None 
INVAL 

This is a dummy subroutine in the main deck which is called by putting a 
minus sign in column 1 of data card 3. (See Data Preparation) It is written 
by the programmer in order to let the computer calculate the initial condi- 


tions of the variables. 
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As an example: 







s + H(21) 


REAL PL 


Figure 2-2(b) Example for INVAL 
The initial value of C(22) has been read in. The sequence for the in- 
itial conditions to be set on varíables C(21) and C(20) is: 


C(21) 
C(20) 


C(22)/G(22) 
C(21)*H(21)/G(21) 


These two cards are inserted in the INVAL subroutine immediately after 
the COMMON cards. 
CURVE 
Usage: CALL CURVE (NO,I,J) NO = 1,2,..,9 
This subroutine allows up to 9 curves, each curve having up to 50 break- 
points. Unlike other subroutines, the first argument, NO, refers to the 
number of the curve vice the box number. For example; 

CALL CURVE (2,7,22) 
generates a value for C(7) as a function of variable C(22) according to curve 
2. The set of breakpoints, x], У, х), Уд» хр, Ук are entered as data 
(See data preparation) and a continuous curve is generatedby a series of 


straight line segments. 
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CURVE 2 


A eorum 


Figure 2-3 Example for Subroutine CURVE 
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2.4 Data Preparation 


The input data required is as follows: 


1. 


2a. 


2b. 


D 


E 


FORMAT (I3) The number of block data cards. One card with 

a O0 or - in colum 1 followed by a two digit number in col- 
umns 2 and 3. The minus sign indicates that there is curve 
data to be read in. 

FORMAT (I2, 5E10.4) Block data cards in any order. Each card 
contains the block number in columns 1 and 2 followed by the 


input data, G(N), H(N), P(N), Q(N), R(N), for that block. 


omit if no curves 
FORMAT (Il) The number of curves. 
FORMAT (I2) The number of breakpoints for Curve l. 
FORMAT (2E20.5) The coordinates of the breakpoints for curve 
1 (Χι, Yi» 7223132) Ух). One card for each coordinate. 
The number of cards must equal the number of breakpoints. 
Continue data cards 2b and 2c for each curve. 


FORMAT (13) Number of initial condition cards. One card with 
a 0 or a - PR 1 followed by a two digit number in 
columns 2 and 3. The minus sign will call INVAL, a subroutine 
written by the programmer which automatically sets the initial 
conditions of the variables. The use of INVAL is optional. 
FORMAT (12,Е20.6) Initial condition cards in any order. 
Each card contains the variable number in columns 1 and 2 
followed by the initial condition for that variable. 


FORMAT (F10.5) Time increment for computation. 
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FORMAT (F10.5) Total time of run. 

FORMAT (F10.5) Time interval between periodic print outs. 
FORMAT (F10.5) Time interval between periodic graph 
points. A maximum of 900 graph points will be plotted. 
FORMAT (212) Variables to be printed. First variable 
number to be printed in columns l and 2. Last variable 
to be printed in columns 3 and 4. The first, last, and 
all inclusive variables will be printed according to the 
time interval between periodic printouts. 

FORMAT (I2, 6A8) Number of graphs desired in column 2 (5 
maximum). Name, etc. in the next 48 columns. This informa- 
tion will appear on every graph. 

FORMAT (212,6A8) One data card for each graph desired. 

A two digit number in columns 1 and 2 for the variable 
which is to appear as the X coordinated, and a two digit 
number in columns 3 and 4 for the variable which is to 
appear as the Y coordinate. ТЕ time is desired as a co- 
ordinate 99 is to be used as the variable number. Graph 
title, etc., in the next 48 columns. This information 


appears only on the graph indicated on that data card. 
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2.5 Example 

In order to show the validity of such a digital computer program a 
third order system with two nonlinearities was chosen and set up on the ana- 
log computer [| 5_Jand on the digital computer using program Analog. The 
system chosen was critical in that the response was unstable for a voltage 
input of 10.2 volts. The amplifiers in the acceleration and tachometer 
feedback channels were set to saturate at three volts. The three volts 
on the digital was an exact quantity, whereas the three volts on the 
analog was nominal as one would expect using diodes. The Analog digital 
program block diagram is shown in figure 2-4 and the analog computer cirucit 
diagram is shown in figure 2-5. Subroutine Diagram and the data for pro- 
gram Analog is shown in figure 2-6. The comments in figure 2-6 explain the 
meanings of the data cards and are not part of the data. Figures 2-7 and 
2-8 are the time responses of the Analog program and the analog computer 
respectively. The two traces are almost identical with the exception of 
the output of the saturated amplifiers, and the reason for that is ex- 


plained above. 


19 





πϑη5ά5 ләрло ратиз эташехя 4-2 AMBTJ 


a era —— - -- ——— —— -ө-- - —— A > 





20 





E . 








-- 9 -— 


шәдеде дердо PTU eiduexe лој чотзетпшт$ дезпашоо Zo[euy $-2 eJn3rJj 


"m 


poo 3 


GBR, 
ya 
F g- 1) "X T 

M | | / 


> 
Го 
х Аа рәзеотрит Er 
әсе э3е3Тол uorygemges" - 


рәлтѕәр доу 395 5304 





Speadejodoru ит O Т 
Suyo3au ит Ц. 
:930N = 
5% + 
POMA y r o > 
τ | | % T к We 
y 
7 7 


1 
< 
+ 
п 
yr: 
<> 


< 
Юю ΠΠ αν 
өң “55 A | + 
I 


21 





'эташехе uəqsÁs JI9pJo pJrU4 лој езер pue WYunVIT euranoadqng 9-2 SMMITA 


В i-r cb 2066 


SHAVED BWIL SA (iNdino) VLBHL 9066 
30 πι ОЛ ео ГЕ 8066 
. C us ВИТЬ сл а | I 6066 
ST SA 100 152 51811 207 2552 
104 аты снаумӘ 3Л1- ЕОС Е ае а ЕИ Е ое a = 
(23S1*) 8 Ol 2 S3'18VINVA іміне 802 
335 £Z0* AH3A3 S3J1HVISVA LONE 20° 
2381* A33A3 SITOVIVVA 3Η1 iNIäd : l° 
ΝΩΝ JO H19N31 vo 
s INIR ONI ἕως” 
SNCILIGNOD IVILINI ON | 
οι SI aA Tv ATS οι 047 8 
6 ONY 2 $32078 €- 1111} 83M07 eç- sero O*l 6 
€ L1IWI] H3ddfl O*IsNIVO eç- "C ele 2 
931N1 9343 ο ο O ST 9 
*1-(650)H *T=(50)9 moe *1 G 
“2-("СУН %09-(%0)9 a *09 5 
/ ОТ 30 43165 ШО! 1 
VIVO SLIM 5ΧΠΟΠΕ JO SIENNN L 
S 1N3NWO2 А 
QN3 
: QN3 
КЕПТЕР 


Ва 
(8g*6*6)H3lIWI' | TIVO 
(5:98*8)VAIS3Q 11V) 

τ: η, ΘΟ 

Wr en oa a 11.2 
узи У 
о DEO ην 
.0.6 1,68 - SATZ 1 
Селам. 15 

.va5viQ 3NILNOYENS 


22 




















E L 


Figure 2-7 Program Analog results for third 
order system example 
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2.6 Conclusion 


Examination of the response curves of program Analog indicates that 
the digital computer program accurately simulates an analog representa- 


tion of a control system. 
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3. The Bailey Combustion Control System. 
3.1 Discussion 

The Bailey Combustion Control System is utilized on the DLG-9 boiler 
which this thesis investigates and analyzes. This section describes how 
the combustion control system operates, and gives general descriptions of 
various components of the system, which will be further amplified in 
greater detail in later sections. 

The boiler steam drum pressure is measured in the steam pressure 
transmitter that develops an air loading pressure proportional to steam 
pressure. Since the pressure transmitter is measuring steam pressure 
over a range from 900 psig to 1500 psig, it will produce an air loading 
pressure over a range of three to 27 psig in direct proportion. Since 
loading pressure is proportional to steam drum pressure it can be substitut- 
ed for steam pressure, and become the steam drum pressure measurement in 
the control system. 

The steam flow is measured by a flow transmitter which develops an 
air loading pressure that is directly proportional to steam flow. Since 
the DLG-9 boiler has a full power generating capacity of 166,000 pounds 
of steam per hour, the flow transmitter will vary its output loading 
pressure from three to 27 psig as the steam flow is changed from zero to 
166,000 1bs/hr. If the boiler is generating 50 percent of full power 
(83,000 1bs/hr), the transmitter will send an air loading pressure of 15 
psig (50 percent scale). 

The amounts of fuel, air, and feedwater required depend upon two facts: 
(1) the amount of steam being used and (2) the desired pressure to be 
maintained. This designates steam flow as our index to changes in boiler 
load, and steam pressure as the correction factor to set the desired 
standard. 
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The steam flow signal is transmitted to a master relay where the 
steam flow measurement is exactly reproduced. The steam pressure is trans- 
mitted to the same relay where it is balanced against a reference force 
that represents the desired steam pressure. Any difference between the 
actual steam pressure and the reference set point will have an inverse 
effect on the output of the relay. The output signal from the master re- 
lay becomes the boiler load demand signal calling for a corresponding firing 
rate. 

To illustrate more clearly, suppose the boiler is generating steam at 
25 percent of its rated capacity, and steam pressure is being maintained at 
1200 psig. The steam flow meter is sending nine psig to the master relay, 
which in turn develops the same pressure. Since the steam pressure is 
equal to the desired value, there is no difference, and has no effect on 
the demand signal. The firing rate of the boiler is equal to 25 percent 
called for by the master relay. Now the engine throttles are opened, in- 
creasing the steam flow from 25 to 50 percent of boiler full power. The 
increased steam flow results in an increase in boiler demand signal, this 
is almost instantaneous, however, due to the lag in the response of pumps, 
blowers, and heat transfer the boiler does not immediately increase its 
generating rate. This results in a decrease in steam pressure. The 
difference in actual pressure versus desired reference causes the master 
relay to increase its output a proportional amount which calls for a fir- 
ing rate in excess of the rate called for by steam flow. This is called 
overfiring, and is necessary to overcome the various time delays encounter- 
ed in the control and plant systems. With the firing rate now greater than 
steam flow, steam pressure is more quickly restored to the desired value. 

On a decrease in steam demand, the firing rate is decreased simultan- 
eously with steam flow but the residual heat in the boiler tends to increase 
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the steam pressure slightly. This then calls for firing at rates less 
than steam flow or underfiring. 

It becomes readily apparent here, that deviations in steam pressure, 
during changes in load is nothing more than a function of the characteris- 
tics of the boiler and its auxiliary machinery. 

The master demand signal is transmitted to an air flow demand Standa- 
trol, and through a limiting relay, to an oil flow Standatrol. A Standa- 
trol is Bailey Company's trade name for their standard proportional plus 
reset controller. In the air flow Standatrol the master demand signal is 
balanced against the air flow measurement and the output pressure from the 
Standatrol serves to position the forced draft blower throttles and air 
flow control dampers. The Standatrol is calibrated to maintain a constant 
output pressure as long as the air flow measurement is equal to the demand 
signal. When the demand signal increases calling for more air flow, the 
Standatrol will commence increasing, at a rate depending upon the relay's 
proportional band and reset settings, until the air flow again balances 
the master demand. The air flow is measured by a transmitter taking the 
drop of flue gases across a section of the economizer and extracting the 
square root. This air flow transmitter signal goes to a calibrating relay, 
referred to as the excess air relay, which has remote bias and proportional 
band adjustments. These adjustments permit the operator to regulate the 
amount of excess air for economy haze, smokeless operation, or to blow 
tubes. The output of the excess air relay is transmitted to the air flow 
Standatrol. With this arrangement any master demand value will result in 
the forced draft blowers running at a speed that will provide a corres- 
ponding air flow regardless of the characteristics of the blower throttles 


or the number of blowers in use. 
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The control of the oil flow is similar to the Air Flow Control 
System except for the inclusion of the fuel limiting relay. Since the slow 
inertia of the blowers causes air flow to lag behind the steam flow on an 
increase in load, it would be possible to admit oil to the boiler before 
there is sufficient air to burn it. For this reason the master demand 
signal is piped to the fuel limiting relay along with the air flow measure- 
ment from the excess air relay. The fuel limiting relay is calibrated so 
that its output pressure can never be greater than the air flow measure- 
ment, but on a sudden decrease in load, the output can decrease at the same 
rate as the master demand regardless of how long it takes the blowers to 
slow down and decrease the air flow. The output from the fuel limiting 
relay is transmitted to the oil flow Standatrol whose principle of operation 
is identical to that of the air flow Standatrol. 

On the DLG-9 boiler the fuel flow metering system is of the return flow 
burner type of fuel oil system. In this system, oil is pumped to the burners 
at a constant pressure of 1000 psi and at a relatively constant rate per 
burner. The design of the burner is such, that if permitted, the fuel can 
flow by the atomizer through the return line to the pump. The fuel that is 
not returned is sprayed into the fire box and burned. A control valve is 
installed in the return oil line and by closing or opening the valve, the 
amount of fuel sprayed into the fire box will be increased or decreased. 

The output from the oil flow Standatrol positions the oil flow control 
valve. 

An areameter transmitter is installed in both the supply and return 
fuel lines and each develops an air loading pressure proportional to their 
respective oil flows. These supply and return flow signals are transmitted 
to a totalizing relay where an output pressure is developed that is propor- 
tional to the flow of fuel into the boiler by subtracting the oil returned 
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signal from the oil supplied signal. This fuel burned measurement is trans- 
mitted to the oil flow Standatrol and balanced against the oil demand signal 
from the fuel limiting relay. 

As in the air flow control, a linear relationship between the master 
demand signal and oil flow is maintained regardless of the flow character- 
istics of the burners or the number of burners in use. 

The DLG-9 boiler employs the three-element feedwater control system 
which is the standard metered system that is widely used in the Navy. 

The basic control elements consist of steam flow as the demand index, 

water flow as the response impulse, and boiler drum level as the super- 
visory impulse. The steam flow and water flow transmitters sense the differ- 
entials created by flow nozzles in the steam header and the feedwater header 
and send out corresponding air pressures. These signals are applied to a 
ratio relay which develops a set point signal representative of normal 

steam drum level when the steam flow and water flow are equal.: Any un- 
balance between steam flow and water flow will shift the set point for water 
up or down depending on the direction of unbalance. This signal is applied 
to the feedwater standardizing relay where it is compared to the actual steam 
drum level signal. The steam drum level is measured by a bellows type force 
balance differential transmitter. The transmitter measures level by compar- 
ing the actual drum level against a reference head, and the measurement is 
indicated by a pointer directly powered by the bellows mechanism. The trans- 
mitter also develops an air loading signal that is directly proportional to 
the steamdrum level in inches. 

Any error between the measured drum level and the set point is recog- 
nized in the feedwater standardizing relay which sends a signal corres- 
ponding to this difference to the proportional plus reset feedwater control- 
ler. The signal from the feedwater controller positions the feedwater 
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control valve which admits more or less feedwater to the boiler depending 
on the signal. 

The Bailey Combustion Control System as described in this section is 
in use in the Navy on destroyers and, with some modifications, on Some air- 
craft carriers. As was pointed out in the thesis introduction there re- 
mains much to be done to improve the response of the boiler to sudden 
changes in load especially in the destroyer type ships; this is intimately 
related to the type of combustion control system which controls the boiler. 
Thus by analyzing and investigating the Bailey Combustion Control System 
and the boiler on DLG-9, the authors feel that through this work adaptive 
type control improvements may be forthcoming, and it is a necessary step 


towards the design of a digital computer control system. 
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4. Combustion Air Flow Control System 
4.1 Discussion 
The air flow control system is comprised of a typical feedback system 
including two principal nonlinearities, the steam control valves, and the 
forced draft blowers themselves. The system is composed of: 
a. Proportional Plus Reset Air Flow Controller 
(Bailey Meter Company "Mini - Line" Standatrol) 
b. Proportional Plus Rate Action Controller 
(Bailey Meter Company Rate Booster Relay) 
c. Forced Draft Blower Actuator, Linkage, and Throttle Valve 
(Bailey Meter Company Model AC-44 Control Drive with Mason- 
Neilan Control Valve) 
d. Carrier Corporation Main Forced Draft Blowers 
e. Total Air Flow Transmitter 
(Bailey Meter Company Type CJ-20 Differential Transmitter 
with Type KC-16 Square Root Converter - Transmitter) 
f. Air Flow Calibrating Relay (Bailey Meter Company Model 
AR-40 "Mini - Line" Relay). 
The system is set up as in figure 4-1 with pneumatic signal Pm as the in- 
coming signal, and the quality of air Qa as the output. The quality of 
air is defined as the amount of air flow in cubic feet per minute divided 


by the rated blower capacity (37,800 cubic feet per minute). 
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4.2 Proportional Plus Reset Air Flow Controller 

This Standatrol is a pneumatic relay designed and calibrated to 
receive a loading pressure, Pe, representing the difference between the mas- 
ter demand signal from the steam pressure controller and the signal from 
the air flow measurement. The output of the Standatrol serves to position 
the forced draft blower throttles and the air flow control dampers. 

The proportional plus integral controller is linear and has constants 
which are shown in figure 4-3. Its minimum output is zero psig and its 


maximum output is 30 psig. 
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Figure 4-2 Forced draft blower data 
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4.3 Proportional Plus Rate Action Controller 

This device is intended to provide increased phase margin in the 
blower control loop thus allowing maximum loop gain. The proportional plus 
integral Standatrol is calibrated to maintain a constant output pressure 
as long as the air flow measurement is equal to the demand signal. Should 
the demand signal increase, calling for more air flow, the rate action 
booster relay will commence increasing its output, Pq, at the rate set by 
the proportional band and reset settings until the air flow again balances 
the demand signal. 

Proportional plus rate action devices unfortunately tend to act as 
noise amplifiers, and thus limiting action was encountered in the simula- 
tion; however, by proper choice of the time increment for the digital 
simulation the "noise effects" of the rate action booster relay can be dimin- 


ished, and favorable results obtained. 
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4.4 Forced Draft Blower Actuator, Linkage and Throttle Valve 

The pneumatic signal, Pq, from the rate action booster relay flows 
into the forced draft blower actuator positioning the throttle valve which 
admits the steam, Gb, to the forced draft blower turbine. 

The blower speed versus loading pressure Pq, curve, figure 4-4, was 
obtained experimentally at NBTL [3] and since the limits of the pneumatic 
devices are between three and 27 psig, there are built in limits on the 
blower actuator, resulting in blower speeds not lower than 1800 RPM, nor 
higher than 7000 RPM. 

The actuator dynamics are presumed to be linear for all blower speeds 
but the gain varies as the blower speed, figure 4-2. This non-linearity is 
due to the characteristics of the V - ported blower steam control valve. 

The non-linear curve of blower steam flow, Gb, versus blower speed, N, 
figure 4-5, was obtained from the experimental data of figure 4-2, and 
since blower speed, N, versus loading pressure, Pq, figure 4-4, is avail- 
able, a function generator (curve three of figure 4-3) can be obtained 
that relates blower steam flow, Gb, as a function of loading pressure, Pq, 
figure 4-6. Using this non-linear function generator, the blower steam 
flow feeds through the dynamics of the actuator into the forced draft 


blower. 
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4.5 Main Forced Draft Blower 

The forced draft blowers are non-linear devices due to the driving 
torque and retarding torque being non-linear functions of turbine steam 
flow and blower speed. 

The blower speed, N, versus air flow, Qb, is a non-linear function 
which was obtained from an experimental run on the DLG-9 test boiler at 
NBTL [3] some of the results of that run are tabulated in figure 4-7. 

The sample calculation below illustrates how the blower speed N is 
related to the air flow Qb. 

At 75 per cent of the boiler full power run (run number 421): 

Blower output (air flow of two blowers) - 36,000: CFM 

Boiler steam flow (Gs) - 130,930 1b/hr 

Blower turbine steam flow (2 Gb) - 13,400 1b/hr 


4.51 Determination of Qb by the Definition of the Quality of Air 


36,000 CFM/2 Blowers 


Ls LL =, 


Ов = (2 В1омегѕ) (37,800 СЕМ) .= &67.7% 4.1 
4.52 Determination of Qb by Boiler Steam Flow 


Following the flow of steam through the plant as shown in figures 


4-8 апа 4-3: 

Pgs = (Steam flow Gs)(Steam flow transmitter gain) plus three psi offset. 
-4 

Bes = (130,930 lb/hr)(1.14 x 10 psi/lb/hr) + 3.0 psi = 17.9 55182 


At the steady state value of 75% full power, from figure 4-8: 


Pep = O and Pbq = Pm = Pgs = 17.9 psi 


thus Qa = (Pbq - three psi offset in calibrating relay) 


(Calibrating relay gain)(Air flow transmitter gain) 
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RUN PERCENT [BLOWER BLOWER 
NUMBER I BOILER ¡OUTPUT TURBINE 
STEAM FLOW: 


2Gí lb/hr) 


55366 cruisinr 


01590 55.27 
130930 uo 
167710 ТОТ; 1.6820 
195930 1208 57130 





Figure 4-7 Experimental data from test run on DLG-9 Boiler | 
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Ops iO 


Q, = (1.05)(0.254 psi/%) = 55.8% 4.3 
And from figure 4-3: | 
Q 
¿E 
EN = (2)(.567) = 49.0% 4.4 


The two values of Qb are in close agreement, the latter value of 49.0% 
was used since the air flow measurement was not as accurate as the steam 
flow measurement. 
4.53 Determination of Blower Speed as a Function of Blower Turbine Steam 
Flow 
2 Gb = 13,400 1b/hr 
Gb = 6,700 1b/hr 

Entering figure 4-5 with 6700 1b/hr the corresponding blower speed is 
5100 RPM. 

Using the value of Qb of 49% and the corresponding speed of 5100 RPM, 
function generator (curve ten, figure 4-3) of blower speed N versus air 
flow Qb, figure 4-9, was obtained. 

4.54 Forced Draft Blower Gain 

The transfer function for the forced draft blower as noted in figures 
4-2 and 4-7 are non-linear with the gain and the dynamics being functions 
of blower speed, N, or air flow, Qb, since by figure 4-9 these two para- 
meters are related. 

The gain function generator (curve four of figure 4-3) is a function 
of Qb and was determined from the data given in figure 4-7, as follows: 

Using the data from run number 421: 


Qb = 49.0% 


Gb 6700 lb/hr 


thus blower gain is defined as Qb/Gb, and for this speed 
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49.0% 
Blower Gain = 6700 1b/hr 20810755 4.5 


At the low speeds, since the lowest value of loading pressure, Pq, 
is three psi, the corresponding blower speed is 1800 RPM, which gives a 
corresponding blower turbine steam flow, Gb, of 1440 1b/hr from figure 4-2. 
From experimental test run on the DLG-9 test boiler at NBTL [3] the corres- 
ponding air flow Qb is 6.35%. Thus the function generator relating blower 
gain to air flow, figure 4-10, was obtained. 

4.55 Blower Time Constant 

The time constant function generator (curve five of figure 4-3) is a 
function of blower speed and was determined directly from the data given in 
figure 4-2 in which the time constant © is given as a function of speed, 
thus figure 4-11, 1/4» ав a function of blower speed, М. 

By means of these various function generators the action of the main 
forced draft blowers can be simulated from ten percent to 120 percent of 
boiler full power. 

4.6 Total Air Flow Transmitter 

This instrument is designed to measure the pressure differential of 
the air flow across an orifice or restriction, extracts the square root of 
the differential pressure, and develops a pneumatic loading pressure that 
is proportional to the flow of the air. The air flow transmitter is simu- 
lated as shown in figure 4-8. 

4.7 Air Flow Calibrating Relay 

This relay, operating on the air flow transmitter output signal, Pqa, 
in the feedback circuit of the air flow control loop, is provided to permit 
scaling of th» air flow open-loop gain in order to obtain optimum air to 
fuel ratio. The gain and bias are both remotely adjustable from the engine 


room control console, with the gain having a range from .5 to 20 and the 
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bias having a range of plus or minus three psi. The air flow calibrating 
relay is shown in figure 4-8 with a positive three psi bias applied. 
4.8 Root Locus Studies 

By means of program Root Locus (Appendix II)root locus studies of the 
air flow control system were undertaken to determine the relative stability 
of the system with respect to possible changes in system gain. 

Two studies were made. One was at a blower speed of 3000 RPM, which is 
near the plant's cruising condition; and the other was at a blower speed of 
7000 RPM which is near the 90% full power condition. The gain of the system 
was varied from three tenths of its normal gain to a value of three times 
its normal gain. 

Figure 4-12 is the root locus plot for 3000 RPM. The closed loop 
system is stable for the normal gain and has roots, which are encircled on 
figure 4-12, at: 

-.1205 +]1.22] 

2.1205 -]1:22/] 

-.34&59 +10.04162 
-.3459 -10.04162 
-.8520 

-2.371 

NI O 20 

-3.981 -19.201 

The system is unstable at 1.505 times the normal gain. 

Figure 4-13 is the root locus plot for the air system at /000 ЕРМ. 

The closed loop system is stable for the normal gain and has roots which 
are encircled on figure 4-13 at: 

-0.1555 +j1.243 

-0.1555 -11.243 

-0.7697 

-0.4029 +10 .295 

-0.4029 -10.295 

- 9.929 19.176 


-2.466 
229729 


The system is unstable at 1.699 times the normal gain. 
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Figure 4-12 - Air system root locus with forced 
draft blower speed of 3000 RPM, 
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draft blower speed of 7000 RPM, 
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4.9 Air Flow Control System Simulation 
4.91 System Simulation 

The air flow control system model is shown in block diagram form in 
figures 4-3 and 4-8. The model was simulated using program ANALOG on the 
digital computer. The inputs to the model were steam flow, Gs, and the 
superheater outlet pressure, Po. The inputs were taken from experimental 
data furnished by NBTL [4] . Figure 4-14(a) is the graph of the superheater 
outlet pressure used as an input in the simulation where the steam flow was 
increased from ten per cent to 90 per cent of full power in 23 seconds. 
Figure 4-14(b) is the graph of the superheater outlet pressure used as an 
input in the simulation where the steam flow was decreased from 90 per cent 
to ten per cent of full power in 23 seconds. 
4.92 Results 

Figure 4-15(a) is the actual system response to a ramp change in steam 
flow from the ten per cent to the 90 per cent full power condition in 23 
seconds. The peak air flow of 87 per cent occurs at 29.5 seconds and air 
flow value is 78 per cent at 70 seconds. 

Figure 4-15(b) is the response of the model to a ramp change in steam 
flow from the ten per cent to the 90 per cent full power condition in 23 
seconds. The peak air flow occurs at 28 seconds and has a value of 89.4 
per cent. The air flow is 77.5 per cent at /0 seconds. 

Figure 4-16(a) is the response of the actual system to a ramp change 
in steam flow from 90 per cent to ten per cent full power in 23 seconds. 

Figure 4-16(b) is the response of the model to a ramp change in steam 
flow from the 90 per cent to the ten per cent full power condition in 23 


seconds. 
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4.93 Conclusions 

Figures 4-15 and 4-16 show that the simulation model accurately por- 
trays the air flow control system. 

There is a difference in the starting points of the actual system and 
the model due to the superheater outlet pressures at which the tests were 
started. The model simulations were started with a superheater outlet 
pressure of 1200 psi while the actual system tests were started at a 
slightly different superheater outlet pressures. A change of ten psi in the 
superheater outlet pressure will cause a change of approximately five per 


sent in the quality of the air. 
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5. Fuel Oil Flow Control System 
5.1 Discussion 
The fuel oil flow control system is of the return flow burner type 
of fuel oil system, which is essentially linear in both static and dynamic 
characteristics. The system is composed of: 
a. Proportional Plus Reset Oil Flow Controller (Bailey Meter 
Company  "Mini-Line'" Standatrol) 
b. Return Fuel Oil Flow Control Valve (Fisher Governor Company 
Model YFES-9M) 
c. Supply and Return Fuel Flow Transmitters (Bailey Meter Com- 
pany Model JR-13 Area Meter) 
d. Oil Flow Totalizing Relay (Bailey Meter Company "Mini-Line" 
Relay) 
The block diagram is shown in figure 5-1 with pnuematic input signal, Pdf, 
as the demand index, and the output being the fuel oil flow of the boiler, 
Gf. There is a constant input of 15,900 1bs/hr of fuel oil supplied 
to the system at a constant pressure of 1000 psi and the burners are de- 
signed to return fuel oil to the fuel oil pump and this amount is deter- 
mined by the return fuel oil control valve. 
5.2 Proportional plus Reset Oil Flow Controller 
This Standatrol, like the one in the air flow control system, is a 
typical proportional plus integral controller. The input pneumatic signal, 
Pf, is balanced against a pneumatic reference, Pdf, from the fuel limiting 
selector relay and the Standatrol developes an output air signal, Pdr, which 
is sent as the control signal to the return fuel oil flow control valve. 
5.3 Return Fuel Oil Flow Control Valve 


This V-port fuel oil control valve produces linear return fuel flow 
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with lift, and the variation in firing rate of the burner with changes in 
return fuel flow is linear. This valve is a diaphragm operated valve with- 
out positioner. The valve has a pneumatic signal, Pdr, as its loading 
signal and its output is A Gr, the return fuel flow. 
5.4 Supply and Return Fuel Flow Transmitter 

These area meters measure the supply fuel flow, A Gp, and the return 
fuel flow, A Gr, respectively, by means of a metering valve, and each trans- 
mits a pneumatic signal, Pr from the return fuel flow transmitter and Pp from 
the supply fuel flow transmitter, which is representative of the amount of 
fuel flow. These pneumatic signals are sent to the oil flow relay. 
5.5 Oil Flow Totalizing Relay 

This relay develops an output pneumatic pressure, Pf, that is linear- 
ly proportional to the flow of fuel oil burned. This signal is obtained 
by subtracting the return fuel oil pneumatic signal, Pr, from the supply 
fuel oil pneumatic signal, Pp. This signal, Pf, which represents the 
fuel oil burned, is transmitted to the oil flow Standatrol and is balanced 
against the oil demand signal, Pdf, from the fuel limiting relay. 
5.6 Stability Analysis of the Fuel Oil Control System 
5.61 Discussion of the Problem 

The first computer test runs which were conducted using the transfer 
functions for the fuel oil control system furnished by NBTL [ 1] showed 
that the fuel oil control system was unstable. This, in turn, caused the 
simulation of the complete boiler to be unstable. 
5.62 Root Locus Study of Fuel Oil Control System 

In order to prove that the fuel oil system alone made the complete 
boiler simulation unstable with the data furnished, a root locus of the 
System was derived from figure 5-2. Since 

Gf = Gp - Gr Si 
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СЕ could be fed thru the feedback loop in stead of Gp - Gr. Also, 
from figure 5-2 

Ср = 15,900 - .0445 СЕ 562 
Adding equations 5.1 and 5.2 yields 

1.0445 Gf = 15,900 - Gr 5.3 
Us.ng the above manipulations, figure 5-2 can be reduced to figure 5-3 
in order to produce a characteristic equation for a root locus plot. 
The root locus of the system, figure 5-3, is shown in figure 5-4&. The 
gain of the system is varied from three tenths to three times the gain of 
the system as shown. From this root locus it is seen that the system is 
unstable for the specified gain with a pair of complex roots in the right 
half plane at +.67 + j10.0. ы 
5.63 Supply and Return Oil Transmitter's Dynamics 

The transient response to a step input test which was conducted at 
NBTL [ 1] on the area meter transmitter is shown in figure 5-5. The trans- 
fer function of the transmitter's dynamics was simulated on an analog 
computer and the transient response to a step input is shown in figure 
5-6(a). This response did not compare favorably with the actual test re- 
sponse, figure 5-5. By adjusting the damping coefficient of the quadratic 
I on the analog computer, a more favorable response was attained (figure 
5-6(b)). The second order approximation of the dynamics of the transmitter 
1 


was found to be .03352 + 0.1955 + 1 . A comparison of the responses 


shown in figures 5-5, 5-6(a), and 5-6(b) are tabulated in table 5-1. 
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Figure 5-4 - Fuel oil control system root locus with 
transmitter dynamics furnished by nBTL. 
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Figure 5-6(a) - Area meter transmitter transient response 
to a step input - Analog computer simulation. 
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- Figure 5-6(b) - Area meter transmitter transient response 
to a step input - Analog computer simulation, 
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TABLE 5-1 


Figure 5-5 5-6(a) 5-6(b) 

Dynamic Transfer dc incu c 
Function 

Rise Time (100%) .43 sec .38 sec .46 sec 

Max Overshoot 127 317 13,57 

Time to peak .55 sec .6 sec .6 sec 

Second Crossover ‚73 вес 1.0 sec 1.1 sec 


This table illustrates why this set of dynamics was chosen. 
5.64 Root Locus Results with New Transmitter Dynamics 

Using the new value in the dynamics for the supply and return oil 
transmitters, the root locus of this system was plotted and is shown in 
figure 5-/. Again a range of gains from three tenths of normal gain to 
three times that gain is plotted. For the system at its specified gain 


the closed loop system is stable with the following roots: 


The closed loop system is somewhat oscillatory because of the complex roots 


at -.25 + j9.9, but it is not dominated by these roots. 
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Figure 5-7 - Fuel oil control system root locus with 
transmitter dynamics determined by analog 
computer results. 
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5.7 Fuel Oil Flow Control System Simulation 
5.71 System Simulation 

The fuel oil flow control system is shown in block diagram for pro- 
gram Analog in figure 5-8. This system was tested and simulated in con- 
junction with the air flow control system, since these two are connected 
by the fuel oil limiting selector relay, and the inputs to the combined 
air flow and fuel oil systems are steam flow, Gs, and the superheater out- 
let pressure, Po, as was discussed previously in section 4.91. 
5.72 Results 

The system response to a ramp change of ten to 90 per cent full power 
in steam flow from a test conducted on the DLG-9 test boiler at NBTL [4] 
is shown in figure 5-9(a). For this test the fuel oil flow reached a 
peak of 107% of fuel oil flow at full power (13,000 lbs/hr), and the 
final value was 87%. 

The results from the digital computer simulation shown in figure 
5-9(b) indicate a peak of 95.6% in 30 seconds and a final value of 87%. 
The test data from NBTL is somewhat questionable since the fuel oil flow 
does not follow the air flow response as stipulated by the fuel oil limit- 
ing relay. As was previously discussed, the fuel oil limiting relay trans- 
mits the fuel oil demand signal, which is the lower of two signals from 
the air flow system and the master demand signal. If the fuel oil flow were 
as indicated in figure 5-9(a) that is, leading the air flow, black smoke 
would issue from the boiler due to excessive firing of the burners. 

Figure 5-10(a) illustrates the response of the test boiler at NBTL [4] 
for a ramp change in steam flow from 90 to ten per cent of full power. In 
this test the fuel oil has a minimum value of 18% at 22 seconds and settles 


out to a final value of 25%. 
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The digital computer simulation results, presented in figure 5-10(b) 


has a minimum value of 15% at 22 seconds and has a final value of 14%. 
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6. Boiler Water Level Control System 
6.1 Discussion 
The water level in the boiler steam drum is regulated by a conven- 
tional three-element feedwater control system. The system is designed in 
cascade fashion; that is, the proportional plus reset control mode regulates 
feedwater flow rate to equal steam flow rate while simultaneously maintain- 
ing boiler drum water level at the desired value. This arrangement is pro- 
grammed to produce a smooth transient response which to some degree compen- 
sates for variations in feedwater flow rate caused by the effects of "shrink" 
and "swell" when changing load. The system consists of: 
a. Proportional Plus Reset Drum Water Level Controller 
(Bailey Meter Company "Mini-Line" Standatrol). 
b. Feedwater Control Valve (Bailey Meter Company Model A-FFD-1B) 
c. Water Flow Transmitter and Steam Flow Transmitter. 
(Bailey Meter Company Model CR-166). 
d. Water Flow Feedback Signal Filter (Bailey Meter Company 
Needle Valve Volume Chamber Assembly). 
e. Drum Water Level Transmitter (Bailey Meter Company Model 
1843X). 
f. Steam Flow-Water Flow Differential Relay (Bailey Meter 
Company "Mini-Line" Relay). 
This system is arranged as shown in figure 6-1 with the pneumatic 
signal, Pgs, corresponding to steam flow as the demand index, the pneu- 
matic signal Pl, representing the boiler drum water level as the super- 


visory signal and the water flow, A Gw, as the output of the system. 
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6.2 Proportional Plus Reset Drum Water Level Controller 

This Standatrol differs from the ones in the air flow and the fuel 
flow systems in as much as this is a three element controller in which 
the reference pneumatic signal, Prl, is compared with the pneumatic signal, 
Pbw', from the steam flow - water flow differential relay, and Pl, which 
is sent from the drum water level transmitter, and generates a pneumatic 
error signal, Pel. This pneumatic error signal, Pel, acts as the loading 
pressure for the Standatrol which transmitts a pneumatic signal, Pdw, to 
the feedwater regulating valve. The action of this Standatrol is typical 
of an integral plus proportional controller. 
6.3 Feedwater Control Valve 

This control valve is a typical V-ported positioner equipped diaphragm 
control valve, which has the characteristic underdamped second order fre- 
quency response. 

The pneumatic signal, Pdw, from the water level controller acts as the 
loading signal for the valve and the output is Gw, the feedwater flow. 
6.4 Water Flow Transmitter and Steam Flow Transmitter 

Both of these transmitters are of the same design, and the following 
description of the water flow transmitter applies to the steam flow trans- 
mitter. The pressure differential related to the feedwater flow, Gw, is 
measured across an orifice and the square root of this differential pressure 
is extracted, thus the transmitter develops a pneumatic loading pressure 
Pgw, which is proportional to the feedwater flow. 
6.5 Water Flow Feedback Signal Filter 

This needle valve-volume tank serves the function of an R-C filter in 
the feedback path which is designed to attenuate the feedwater control 
system response, Pgw, to the high frequency components in the measured 


water flow signal, Pgs. 
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6.6 Drum Water Level Transmitter 

The water level transmitter is a bellows type inverse acting differ- 
ential pressure transmitter. It is designed to measure and transmit differ- 
ential head pressure between water level in the boiler drum and a reference 
column of condensate connected into the steam space above the water level 
in the steam drum; the steady state calibration curve changes one psig 
(Pl) for each one inch of boiler level change (A^ L). The pneumatic out- 
put of the water level transmitter, Pl, is sent to the water level controller 
as described in section 6.2. 
6.7 Steam Flow - Water Flow Differential Relay 

This relay subtracts the pneumatic output from the feedwater flow 
transmitter, which has been attenuated through the water flow filter, Pgw' 
from the pneumatic signal, Pgs from the steam flow transmitter and develops 
an output pneumatic pressure Psw' which is linearly proportional to this 
difference in pressure. This pneumatic pressure, Psw', is sent to the water 
level controller as described in section 6.2. 
6.8 Stability Analysis of the Water Flow System 
6.81 Discussion of the Problem 

The first digital computer test runs with the water system isolated, 
using inputs recorded in test runs by NBTL ú indicated that the descrip- 
tion of the system by the transfer functions determined was somewhat oscil- 
latory. These results led to root locus studies of the system. 
6.82 Root Locus Studies of the Water Flow Control System 

The water level and water flow control systems are completely inter- 
connected and the feedback of the water flow system is such that a complete 
root locus study could not be made. Therefore, the effect of the steam 


flow with the water level input held constant was the locus studied, This 
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root locus is shown in figure 6-2. For the parameters of the actual system, 
shown in figure 6-3, the location of the closed roots of the system are: 

ΠΕ 15 

-.042 

es, 

-1.32 + 14.18 

=,405 + 1116 

-29.5 

Even though the pair of complex roots at -.076 + 1.75 are not the dominant 
pair of roots because of the root located at -.042, this pair is highly 
oscillatory and are the probable cause of the oscillation in the feedwater 
responses for the tests conducted on the digital computer. This complex 


pair stems from the feedwater regulating valve and the real root at -.042 


is from the water level controller. 
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Figure 6-2 reedwater control system root locus 
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6.9 Feedwater Flow Control System Simulation 
6.91 System Simulation 

The feedwater flow control system is shown in block diagram for 
Program Analog in figure 6-4. This system was tested and simulated in- 
dependently with data furnished by NBTL ЕЗ . The inputs to the system 
are steam flow and water level as shown in figures 6-5(a) and 6-5(b) as a 
result of maneuvering the boiler from ten to 90 per cent full power in 23 
seconds. 
6.92 Results 

The results as shown in figure 6-6(b) indicate that the system as simu- 
lated is oscillatory for the first 50 seconds and does not agree with the 
test data as shown in figure 6-6(a). The oscillations are due to the root 
locations of the complex pair at -.076 + j./5 as discussed in section 6.82. 
It appears that this transient response is superimposed on the response 
caused by the water level controller and after the initial oscillations 
died out the water flow closely followed the test results. Without re- 
course to the actual test equipment, further study of the parameters of 
this system was not possible and the authors feel that further tests should 


be conducted on this part of the system. 
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7. Computer Simulation of the Complete DLG-9 Steam Generator 
7.1 Discussion 

The complete steam generator, which is composed of the combustion air 
flow control system, the fuel oil flow control system, the feedwater flow 
control system, and the boiler transfer functions is illustrated in block 
diagram form in figure 7-l. This system as shown is in a "delta system" 
form as received from NBTL [1] . The entire steam generator system was 
investigated in the following manner. First the system was simulated 
employing the transfer functions as given for the cruising conditions and a 
small ramp input was used to perturbate the system; next this same pro- 
cedure was used with the 90 per cent full power transfer functions. The 
air flow control system was then changed from the linear condition with 
two operating points to the non-linear combustion air flow system as 
analyzed in section 4 and shown in figures 4-3 and 4-8. The entire plant 
was maneuvered from ten per cent full power to 90 per cent full power using 
first the cruising condition transfer functions for the boiler and then 
using the boiler's 90 per cent full power transfer functions. Finally the 
authors attempted to develop non-linear transfer functions for the boiler 
based on the ones given at cruising conditions and at 90 per cent full 
power conditions as a function of steam flow. 
7.2 Delta System Simulation 

The steam generator system was simulated under two separate conditions 
using the transfer functions given for cruising conditions and 90 per cent 
full power respectively. These were appropriately employed in the forced 
draft blowers and actuators, and in the various boiler transfer functions 
for the steam pressure and water level; the remaining associated equipment 
was simulated as described in the preceding sections. The system was 
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simulated using program Analog which was described in section 2, The sub- 
routine INVAL was the same for both situations and is presented in Appen- 
dix III. This subroutine as described in section 2.3 is used to set up 
the appropriate initial conditions on the various blocks of the system. 
7.21 Delta System at Cruising Conditions 

The subroutine DIAGRAM and the data for the cruising condition simula- 
tion is presented in Appendix III. The steam flow is ramped from cruising 
conditions (56,000 1bs/hr) to 59,000 1bs/hr in five seconds and held at 
this value for 50 seconds as shown in figure 7-2(a). The delta water level 
response to this positive ramp is shown in figure 7-2(b); the peak of .18 
inches is reached in 20 seconds and the level approaches the steady state 
value of zero. The delta steam pressure given in figure 7-3(a) shows a 
dip in steam pressure of a maximum 3.75 psig at 11.5 seconds as expected, 
and the steam pressure approaches a zero steady state value. The air flow, 
figure 7-3(b), exhibits the type of performance one expects for a small 
change in steam flow. The delta fuel flow, given in figure 7-4(a), showed 
some slight oscillations as the delta fuel flow approached its final value 
of 180 lbs/hr. Figure 7-4(b) shows the delta feedwater flow and this system 
as explained in section 6.92 is oscillatory; also the system is inherently 
slow and its final value of 3000 1bs/hr is not reached during the duration 
of this computer run. These results agreed favorably with what was ex- 
pected from actual boiler operation; however there was no test data avail- 
able with which to compare the responses. 

The delta system at cruising conditions was also perturbed with a 
negative ramp in steam flow from cruising conditions (56,000 lbs/hr) to 
53,000 1bs/hr in five seconds and held at 53,000 lbs/hr for 50 seconds as 
shown in figure 7-5(a). The delta drum water level response reaches a 
minimum value of -.18 inches at 12 seconds and then slowly approaches 
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the steady state value of zero, as shown in figure 7-5(b). The change in 
boiler steam pressure, given in figure 7-6(a), reaches a minimum of two 
psig at 36 seconds and then approaches its zero steady state value. Note 
that there is a small rise in pressure at first as expected; but then at 
about ten seconds the pressure goes negative due to the response of the 
feedwater controller, which is somewhat oscillatory (figure 7-7(b)). The 
air flow, figure 7-6(b) exhibits little change as expected for the small 
perturbation of the steam flow. The delta fuel flow approaches its final 
value of -180 lbs/hr with some small oscillations as shown in figure 7-7 
(a). The delta feedwater flow, figure 7-7(b), as mentioned above, was 
oscillatory as in the case of a positive ramp change, and the system is 
slow, not reaching its final value in the 55 seconds of simulation. These 
simulations gave reasonable results with the exception of the feedwater. 
The results, therefore, somewhat proved the validity of the transfer func- 
tions given for the delta plant at cruising conditions for perturbations of 
five per cent around this set point. 
7.22 Delta System Simulation at 90 Percent of Full Power Conditions 

The complete steam generator with associated control system was simu- 
lated at the 90 percent of full power condition using the delta transfer 
functions as given in figure 7-1. The system was perturbed with a five per 
cent ramp in both the positive ànd negative directions as in the previous 
section. The positíve ramp in steam flow from 152,000 1bs/hr to’160,000 
lbs/hr in five seconds and held at 160,000 1bs/hr for 50 seconds is shown 
in figure 7-8(a). The delta drum water level hits a peak of .3 inches in 
13 seconds and then approaches the steady state value of zero as shown in 
figure 7-8(b); the swell observed here is the result of the sudden evapora- 


tion in the riser due to the reduced drum pressure caused by the sudden 
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increases in steam flow. This phenomenon is also observed in the delta 


system at cruising conditions when the positive ramp change in steam flow 
is applied. Figure 7-9(b) shows that the air reaches its final value of 

84 per cent quality in 20 seconds and remains at that value. The change in 
fuel flow, figure 7-10(a), has a slight overshoot and arrives at its final 
value of 640 lbs/hr in 25 seconds. The delta feedwater flow shown in 
figure 7-10(b) still exhibits the oscillatory nature due tv the feedwster 
level control system as explained by the root locus study in -ection 6.92, 
and the change in feedwater is approaching its final value of 8000 Ibs/hr 
and would probably steady out at this value about two minute: softer the 
change in load. 

A negative ramp in steam flow was applied to the delta system at 90 
per cent condition which lowered the steam flow from 152,000 Lb=/hr to 
144,000 lbs/hr in five seconds; this steam flow is shown in figure 7-ll(a). 
The change in water level to this negative ramp is given in figure 7-11(b); 
the drum water level reaches a minimum of -.38 inches in ten seconds, then 
approaches a zero steady state value. The response of the delta steam 
pressure was oscillatory and negative, as shown in figure 7-12(3); this is 
due in part to the oscillatory nature of the feedwater control system as 
explained in section 6.92, The steam pressure should have had a small 
positive overshoot in the first ten seconds and then crossed the zero axis 
and had a negative delta pressure as it approached its zero final value. 
The air flow, given in figure 7-12(b), exhibits the expected response settl- 
ing out at a final value of about 75% in 30 seconds. The change in fuel 
flow is given in figure 7-13(a), which shows that the delta fuel ficw settles 
out at -750 1bs/hr in about 35 seconds. The response of the delta feed- 


water flow was again oscillatory as it had been in all cther rurs, thia 
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simulation was for a period of 55 seconds; therefore, the deita feedwater 
flow was approaching, but did not reach its final value of -8000 1bs/hr 
in this time as shown in figure 7-13(b). The delta system at 90 per cent 
conditions with the negative ramp exhibited the expected transient res- 
ponses with the exception of the delta steam pressure response as ех- 
plained above. 


7.3 Simulation of the DLG-9 Steam Generator with the Non-Linear Air Flow 
System 


The steam generator system was simulated, with the non-linear air flow 
control system shown in figure 4-3, as explained in section 4.91, replac- 
ing the linear system at two operating points as in the delta system simu- 
lation of the previous section. The steam flow was ramped from ten to 90 
per cent of full power in 23 seconds as shown in figure /-i4(b), which is 
the same as the steam flow ramp used by NBTL [4 | in conducting their test 
runs on actual boilers shown in figure 7-14(a), and the responses from 
these test runs conducted by NBTL were used as the reference for compari- 
son with the digital computer simulation results. The digital computer 
simulations were run using three different transfer functions for the boiler, 
namely the linear delta transfer function for the crusing conditions, the 
linear delta transfer function for the 90 per cent full power conditions, 
and lastly non-linear delta transfer functions for the boiler which are a 
function of the steam flow. These last transfer functions for the boiler 
are generated by using the cruising conditions transfer functions for values 
of steam flow below cruising conditions; and for values between cruising 
conditions and 90 per cent full power the gains and time constants of the 
boiler transfer functions are given a value as a function of steam flow us- 
ing a straight line approximation between the values at cruising conditions 


and 90 per cent full power. 
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The air flow comparisons are shown in figures 7-15(a) through 7-15 
(d), where figure 7-15(a) is the air flow test data from a run on the 
DLG-9 test boiler at NBTL [4] . The air flow figure 7-15(b), obtained 
using the cruising conditions transfer functions for the boiler, hits a 
peak of 80% at 25 seconds as compared with 877 in 30 seconds from the test 
data; also the air flow from the test data settles out at 787 of full power 
at about 45 seconds, whereas the air flow of figure 7-15(b) has a value of 
70% at 50 seconds. The air flow obtained using the 90% transfer functions 
for the boiler, shown in figure 7-15(c), has no overshoot but rather it 
reaches its final value of 90% of full power at about 30 seconds. This 
final value is higher than the test data results of 78%. The air flow 
obtained by using boiler transfer functions which are a function of the steam 
flow is given in figure 7-15(d). The response has no overshoot and settles 
out at about 90% of full power in 28 seconds. This result is of the same 
shape as the air flow using the 90% of full power boiler transfer functions. 
It should be noted in the figures that the air for the test boiler has an 
initial value of about 20%, whereas the air for the simulations on the digi- 
tal computer has an initial value of 10%. 

The fuel oil comparisons are shown in figures 7-16(a) through 7-16(d). 
The fuel oil flow from NBTL test boiler for the 10-90% of full power run is 
shown in figure 7-16(a). It hits a peak of 108% at 30 Seconda and settles 
out to 88% of full power in 70 seconds. Using the cruising conditions 
transfer functions for the boiler in the simulated run, the fuel oil hits 
a peak of 88% full power at 25 seconds and settles out to a final value of 
70 per cent in 70 seconds. Using the 90% boiler transfer functions the 
fuel oil flow hits a peak of 102% at 39 seconds and has a value of 90% 
of full power at 50 seconds. This response is shown in figure 7-16(c). 
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The response of the fuel oil flow using the boiler transfer functions as 
a function of steam flow is shown in figure 7-16(d). This response has no 
overshoot and settles out to a value of 102% after 40 seconds. 

The water flow responses are shown in figures 7-17(a) through 7-17(d). 
The test run response of the actual boiler is shown in figure 7-17(a). 

This response shows that the system is extremely slow with no oscillations. 
The responses from the simulated boiler show that there is one or several 
transfer functions for this system in error or that the method of simulation 
is in error. In all three cases the water response neglecting the oscilla- 
tions is somewhat close to the actual, Using the cruising conditions boiler 
transfer function, the water response, shown in figure 7-17(b), has a value 
of 36% at 50 seconds, whereas the actual response of the boiler, figure 7-17 
(a), has a value of 60%. With the 90% conditions after 50 seconds the water 
flow has a value of 56% as shown in figure 7-16(c); and with the boiler trans- 
fer functions as a function of steam flow, the water response, figure 7-17 
(d), has a value of 50% at 50 seconds. 

The superheater outlet pressure for the actual boiler is shown in figure 
7-18(a), and the outlet rescue for the .simulated boiler using the transfer 
functions as a function of steam flow, is shown in figure 7-18(b). The 
actual boiler has a minimum pressure of 1115 psig, whereas the simulated 
boiler has a minimum of 349 psig. Using the cruising conditions transfer 
functions, the minimum pressure is 1076 psig; and using the 90% full power 
transfer functions the minimum pressure is 1081 psig. These minimums for 
the actual and the simulated boiler all occur at 23 seconds. 

The water level of the actual boiler, shown in figure 7-19(a), has a 
peak of 6.0 inches. For the simulated boiler using the cruising conditions 
transfer functions, the peak water level was 7.7 inches; and for the 907 
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full power conditions the peak was 4.8 inches. These peaks of water level 
for the actual and for the three various Simulatoren occurred at 25 
seconds. 

7.4 Discussion of Results 

The test results for the simulation of the complete boiler with the 
boiler transfer functions dependent upon the steam flow of the boiler .show 
that this is not a good simulation of the actual boiler. The pressure changes 
are excessive and other responses in general do not compare favorably to the 
responses of the actual boiler. 

Simulation of the boiler using the cruising conditions for the boiler 
transfer functions may be accurate in the range of five per cent around 
the cruising conditions set points, but is not very good for large pertur- 
bations, namely from ten to 90% full power as noted by the results of these 
simulations. 

For large changes in steam flow the 90% full power transfer functions 
appear, from the results, to be the best data available with which to simu- 
late the complete boiler for large er О 
7.5 Recommendations 

The authors recommend that the original data recorded for the dynamics 
of the water E at NBTL [1] and the simulation attempted in this thesis 
be closely compared in order to find the cause of the oscillations in this 
simulation. 

It is further recommended that data be obtained at ten per cent full 
power and at about 60 per cent full power in order to have a more accurate 
description of the boiler in a change from ten to 907 full power. 

From the digital computer program results,using boiler transfer 


functions as a function of the steam flow of the boiler, it was noted that 
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the steam pressure/steam flow section of the boiler was the major cause 


of the low superheater outlet pressure developed. 
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Appendix I 


Program ANALOG 


ll; 





о М 


N Y > 


22. 


700 


701 
200 


209 
1205 
205 
204 
203 
202 
201 
9001 
900 
901 


CALL ANALOG 
END 


ine ANALOG 
DIMENSION 6 (99) »Н( 99) ›Р (99 } ›С@ (99) ›В (99) ›С (99 ) › СО (99) › СОО (99) »› 
1С000(991›С0000(99) ›Т (20) »Х1(900) »Х2 (900) »Х3 (900) »Х4 (900) »›Х5 {900 } » 
2Ү1 (900) »Ү2 (900) ›ҮЗ (900 ) »Ү4 (900) »Ү5 (900) 2 (26) »_ (26) » 
3ITITLE(1g)» JTITLEC 12) sKTITLEC12) s LTITLE£ 12) $ MTITLE(12)» 
4 Х(9,50)% Ү(9»50) 
СОММОМ бэНэ»Р»0»ЕБ»С»С0»С00»С000»С0000»97%Х19Х29Х3»Х4»Х5»9Ү1»Ү2»ҮЗ» 
1 Ү8»Ү5,2»)| »Х»Ү 
DO 1 M=1999 
С(М)-0.0 
СО(М)-:0.0 
СОО(М)-0.40 
.СО00(М)-0.0 
С0О000(М)-0.0 
DO 2 N=1»20 
T(N)20,0 
С(99)-:0.0 
LLzÓ 
DO 22 КК-700»900 
Х5(КК)-0.0 
У5 (КК) =О. 
CALL INITCON 
READ IN THE VARIABLES TO BE PRINTED 
FORMAT (I2) 
READ 7OOsNVAR 
KNVAR=NVAR+1 
DO 701 K=2»KNVAR 
READ 7005L(K) 
READ IN THE NUMBER OF GRAPHS 
READ IN THE GRAPH DATA AND TITLES 
FORMAT (1296A8) 
READ 200sNGRAPHSs (I TITLE( I) 51.6) 
IF (NGRAPHS) 9001590015209 
GO TO (201>202>203»204+205) »NGRAPHS 
ЕОРМАТ(212»6А8) - 
READ 1205»К5915»(МТІТ(Е(1)»1-7»12) 
READ 1205sK4sL4e(LTITLE( 1) 9 1=7912) 
READ 1205oK3eL3e(KTITLE(1) »›1=7›12) 
READ 1205sK2sL2s(JTITLE(I) 9127912) 
READ 1205 »K1»L1s(ITITLE(I)»1I=7»12) 
PRINT 900 
FORMAT (6X s4HT IME» 30X sSHVARIABLES ) 
FORMAT (200X951 20) | 
PRINT 901s(LIN)»N=2»KNVAR) 
STORING POINTS TO BE PLOTTED 





299 
210 
211 
812 
810 
811 
504 


503 
502 
501 
500 


800 


799 
705 


735 
734 
733 
732 
731 
που 
729 
728 
mad 
726 
725 
724 
123 
722 
721 
720 
719 
718 
717 
716 
715 
714 
713 
τις 
711 
902 


903 
904 


e — — ed 


C(99) IS PROBLEM TIME 

CALL DIAGRAM 

IF(NGRAPHS) 80058005211 
IF(C(99)-T(6)) 80058125812 
IF(901-LL) 80058005810 

LL=LL+1 

GO T0(5009501 950295039504) sNGRAPHS 

X5(LL)=C(K5) 

Y5(LL)=C(L5) 

X4(LL)=C(K4) 

Y4(LL)=C(L4) 

X3(LL)=C(K3) 

Y3(LL)=C(L3) 

X2(LL)=C(K2) 

Y2(LL)=C(L2) 

X1(LL)=C(K1) 

Y1(LL)=C(L1) 

T(6)=T(6)+T(4) 

CONTINUE 

PRINTING ROUTINE 

IF (NVAR) 73697369799 

IF (C(99)-T(8)) 704s 7059705 

GO T0(7119712971397149715971697179718971997209721972297239724>5 


1725972697275 ар 731,732, 733, 734,735) »NVAR 


2(26)=C1L(26)) 
Z(25)=C(L(25)) 
2024) =Сс(1 (24) ) 
2023) =С(1 (23) ) 
Z(22)=C(L(22)) 
Z(21)=C(L(21)) 
2(20)=C(L(20)) 
Z(19)=C(L(19)) 
Z(18)=C(L(18)) 
Z(17)=C(L(17)) 
Z2(16)=C(L(16)) 
Z(15)=C1L(15)) 
Z(14)=C(L(14)) 
2(13)=C(L(13)) 
2(12)-С(((12)) | ; 
Z(11)=C(L(11)) " 
Z(10)=C(L(10)) 
Z(09)=C(L(09)) 2 
Z(08)=C(L(08)) ` 
Z(07)=C(L(07)) 
Z(06)=C(L(06)) 
Z(05)=C(L(05)) 
Z(04)=C(L(04)) 
2(03)=C(L(03)) 
Z(02)=C(L(02)) 
FORMAT(F15.9) 
PRINT 902sC (99) 
FORMAT (20X»5E20.6) 
PRINT 903, (2 (3) №2, КММАК) 
FORMAT (///) 
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910 


704 
736 
2000 


1000 
850 
851 
865 

1865 

2865 
864 

1864 

2864 
863 

1863 

2863 


1862 
2862 


861 
890 


G 
10 
20 
21 
τ 
50 
51 
60 
c 


4 0 
45 


PRINT 904 

T(8)=T(8)+T(3) 

CONTINUE 

IF(T(2)-C(99)) 10005,100052000 

С(99)-С(99)4Т(1) 

60 ТО 299 

( АВ-ан 

IF(NGRAPHS) 8905,890,851 

GO T0O(861,8625,863,864,865) ,NGRAPHS 

DO 1865 1=1,6 

MTITLE(I)=ITITLE(I) 

CALL DRAW(LL ›Х5 ›\5 ›О ›0 › АВ »MTITLE»0,0»,1»1»0»0,7»2 sl LAST) 
DO 1864 :1=1,6 к 

LTITLE(II=ITITLE(ND 

CALL DRAW(LL 9X45Y490905LAB sLTITLE»0»,0»,1»1»0,0, 7,2 »1»LAST) 
DO 1863 1=1,6 

KTITLECI)=ITITLECI) 

CALL DRAW(LL»X35,Y350505LAB »KTITLEsQ5051»15050»752 sl LAST) 
DO 1862 Iz156 | | 
JTITLE(I)=ITITLE(I) 

CALL DRAW(LL»X2,Y2»0,0»LAB 9ТІТІЕ»0»091>1»0»0»7»2 »1»LAST) 

MOD=0 

CALL DRAW(LL»X1»Y1»MOD,O,LAB sITITLEs090919120909798 »1»LAST) 
GO TO 4 

RETURN 

END | 

SUBROUTINE INITCON 

DIMENSION 6(99)»›Н( 99 ) »Р (99 } ›О (99 } ›В (99) ›С( 99 } ›СО( 99) ›СОО (99) › 
1С000 (99) »С0000 (99) »›Т( 20) »Х1 (900) ›Х2 (900) ›ХЗ (900) »Х4( 900 } »Х5 (900 ) › 
21 (900) »\2( 900) »УЗ (900 } »›\4 (900) ›\5 (900) ›2(26)»[(26)» 
ZITITLE(12) sUTITLE(12) »5KTITLECI2) LTITLE(C12) »MTITLE(C12)» 
4 Х(9»50)»Ү(9»50) | 

СОММОМ бэН»Р»0»К»С>СО»С00»С000»С0000»9Т7»9Х1»9Х2»Х3»Х4» Х5»Ү1»Ү2»>ҮЗ3» 
] Y4»Y55»Z»UL»X»Y 

READ THE NUMBER OF BLOCKS AND CHECK FOR FUNCTION GENERATOR 
FORMAT (1H1) 

PRINT 10 

РОВМАТ ( 13) 

READ 20»NB 

NA=XABSF (NB) 

FORMAT(2H N»10X»2H G6»13X»2H H»13X5»2H P»13X»2H Q»13X»2H R) 
PRINT 21 

READ IN THE BLOCK DATA 

FORMAT(1295E10¢4) - 

DO 60 Jz1»NA 

READ 50sN» IGUN) HUN) SPUN) SQUN) SRCND ) 

FORMAT(1H0,12»5E15.5) , 

PRINT 51, №» (№) НЕМ) эР (№ за (М) RIND 

PRINT 10 

TEST TO SEE IF NB WAS MINUS 

IF (NB=NA) 30540930 

PRINT 45 

FORMAT( 34H THERE ARE NO CURVES TO BE READ IN) 

GO TO 100 
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30 
70 


81 


85 
86 
87 
80 
100 
110 
111 


115 
120 


145 
140 


130 
149 
147 
150 


116 
155 


160 


СУ СУ СУ 


170 


200 


READ IN THE NUMBER OF FUNCTION GENERATORS 
КЕАО 70»МҒС 
FORMAT (Il) 
READ THE DATA FOR EACH FUNCTION GENERATOR 
DO 80 K=1,»NFG 
READ IN THE NUMBER OF BREAKPOINTS 

PRINT 10 
FORMAT (12) 
READ 81 ,NBKP 
READ IN THE BREAK POINTS 
FORMAT(2E20.5) 
READ 85 (X(KeJ) sY¥(KeJ) sJ=1 »sNBKP) 
FORMAT i(13H CURVE NUMBER» I?) 
PRINT 86К | 
FORMAT(I12sX»12HBREAK POINTS) 
PRINT 87, NBKP 
РВІМТ 85 (Х(К»./)»Ү(К»./)».)з1»МВКР) 
READ THE NUMBER OF INITIAL CONDITION CARDS 
PRINT 10 
READ 110»NINC 
ҒОВМАТ(12) 

IF (NINC) 11151405111 
KNINC=XABSF (NINC) 
READ IN THE INITIAL CONDITIONS AND PRINT THEM 
DO 120 M=1»KNINC 
ҒОВМАТ(12%»Е20.6) 
READ 115,N»C(N) 
CHECK TO SEE IF INVAL ROUTINE IS TO BE CALLED 
IF (NINC-KNINC) 130914059130 
FORMAT(38H NO ENTUM. VALUE SUBROUTINE WAS USED »////) 
PRINT 145 
GO TO 149 
CALL INVAL Ä 

PRINT 147 
FORMAT(28H CIN) INITIAL CONDITIONS ) 

DO 160 N=1,99 , 

ІЕ(С(М))155»9160»9155 
ҒОВМАТ(2Х%12>Е20.6»//) 
PRINT 116+N»C(N) 
CO(N)=C(N) 6 i 
COO(N)=C(N) 
COOO(N)=C(N) 1 
COOOO(N)=C(N) | 
CONTINUE 
READ IN THE TIME INCREMENTS T(1) IS THE TIME INCREMENT,T(2) IS TH 
TOTAL TIME OF THE RUN»T(3) IS THE TIME INCREMENTS FOR PRINT OUTS, 
T(4) IS THE TIME BETWEEN PLOT POINTS 
FORMAT(F10.5) 
ВЕАО 170»Т(1) 
READ 1709T(2) 
READ 1709T(3) 

READ 1709T(4) 

ҒОВМАТ(/////) 

PRINT 200 
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210 FORMAT (9X99H TIME INC»11X»14H LENGTH OF RUN»6X» 
1 15H PRINT TIME INC»5X»14H PLOT TIME INC ) 
211 FORMAT (4Е20.8) 
212 FORMAT(//) 
PRINT 210 
PRINT 212 
PRINT 211» T(1),T(2),T(3),T(4) 
PRINT 10 
RETURN 
END ` 
SUBROUTINE INVAL 

ОІМЕМ5ТОМ 6(99)»Н(99)»Р(99)»0(99),К(99)»С(99)»С0(99) »С00(99)» 
1С000(99) »С0000(99) »Т(20) »Х1(900) »Х2(900) Х3(900) »Х4(900)»Х5(900) 
2Y1(900)5»Y2(900)5»Y3(900)»Y4(900),Y5(900) ,2(26) ,L(26)» 

ZITITLE(12) sUTITLE(12) sKTITLE(12) sLTITLE(12) 3 MTITLE(12)» 
6 Х(9,50)»Ү(9»50) 

СОММОМ 6»Н»Р,0О»В»С»С0»С00»С000»С0000>Т» Х1»Х2»Х3%Х4»Х5»Ү1»Ү2»Ү3» 

1 Ү4»Ү5»224»9Х>Ү 
Р(99)-0.0 * 
RETURN 
END 

SUBROUTINE ADVANCE (N I J) 

ОІМЕМ5ІОМ 6(99)»Н(99)»Р(99)»0(99)»Е(99) С(99)»С0(99)»С00(99)» 
1С000 (99) »СОО00 ( 99} »Т( 20) »Х1 (900) ›Х2 (900) »ХЗ (900) эХ4 (900) »Х5 (900 ) 
2\1 (900) »\2 (900) ›\З (900 ) »\4 (900) »\5 (900) »2( 26)» (26) о 
31TITLE(12)»3TITLE(12)»KTITLE(12)»LTITLE(12)9»MTITLE(12)» 
 Х(9,50)»Ү(9»50) 

СОММОМ С»НЭР» vi Е,С>СО»С00» СО00» С0О00О» 7»Х1»Х2»Х3»Х8»Х5»Ү1»Ү2»Ү3» 
1 Ү4,Ү5,794»Х 

AT 

COOO(I)sCOO!I) - 

COOlI)sCOlI) ^ 

CO(1)=C(I1) 

C(I)=1, 4375%С (271%. 25XC0lJ)-«625*C00(J)-425*CO00(J)4. 1875*COOOOt J 

RETURN 

END 

SUBROUTINE DERIVA(Ns I 9J) 

DIMENSION 6(99)»Н(99)»Р(99)»0(99),К(99),С(99),С0(99) С00(99)» 
1С000(99)»С0000(99)»Т7 (20) Х1(900)»Х2(900) Х3(900)»Х4(900)»Х5(900) 
2\1 (900) »У2(900) ›УЗ (900 } »\4 (900) »\5 (900) »Ζ(26!) »ἱ (26! » 
317TITLE(12)>»JTITLE(12)>»KTITLE(12)»LTITLE(12)9MTITLE(12)» 

4 Х(9»50)»Ү(9»50) 

COMMON бэН»Р»0»В»С»С0»С00»С000»С0000»Т»9Х1»Х2%9Х39Х4»Х5»9 Ү1»Ү2»Ү3» 
1 YasY5sZsLsXsY 

С0000(1)-С000(1) 

СООО(1)=СОО(1) 

СОО(Т) «СО(Т) 

CO(I)=C(I) 

C(1)=(G(N)/T(1))*(+5*CO0(J)-— 2.0%С0(.)%1.5%С(,))) 

КЕТОЕМ 

ЕМО 

SUBROUTINE INPUT(N» I] 

DIMENSION 6(99),H(99) )P(99)»Q( 99) 9R(99) ,C(99) ,CO( 99) COO( 99), 
1€000(99) »C0000(99) 9 T(20) 9X1(900) 9X2(900) 9X3( 900) 9X4( 900) 9X5(900) 
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2\1 (900) »\2 {900} »\З (900 } э\4 (900) ›\5 {900} ,2 ( 
3ITITEECI2) s JTITLEC I2  sKTITLEUI2) SL TITLE( ( 12 
4 X(9550)»Y(95,50) 

COMMON ОэН»Р»0»К»С»СО»С00»С000»С0000»7»Х1»Х2»Х3»Х4»Х5»Ү1>Ү2»Ү3, 
1 Ү84»Ү52>24»ХҮ 

COOOO(I)2COOO(lI) 

СООО(І)-СОО(1) 

СОО(1,)=СО(1) 

CO(1)=C(1) 

G(N) IS THE VALUE OF THE STEP,H(N)-THE SLOPE OF THE RANP» 

P(N)-AMPLITUDE OF THE SIN WAVE sQ(N)—-W( FREQUENCY ) OF THE SIN WAVE 

IF (QUN)) 19291 

SINE=P(N)*SINF(QLN)*C(99)) 

GO TO 3 

SINE=0.0 

C(I)=GUN)+H(N)%*C(99)+SINE 

RETURN 

END 

SUBROUTINE REALPL(N»I»J) 

DIMENSION С(99)›Н(99)оР (99) $30 (99) ›К (99) »С (99) ›СО( 99) ›СОО ( 99} » 
1С000(99) С0000(99)»Т(20)»Х1(900)»Х2(900)»Х3(900)»Х4(900)»Х5(900) 
2\1 (900) »\2 ( 900} э\З (900 } »\4 (900 } »\5 (900) ›2(26)› (26) 
31TITLE(12)»JTITLE(12)»KTITLE(1219LTITLE(12)»MTITLE(12)» 

4 Х(9,50)»>Ү(9»50) 

COMMON б»Н»Р»0» ЕВ» С»СО»С00»С000»С0000»7%9Х1»Х2»Х3»Х4»Х5%Ү1»Ү2»Ү3; 
1 Ү4»Ү5,2%9(»Х>Ү 

С0000(1)-С000(1) 

COOOlI)sCOO(I) . 

COO(I)sCOL(lI) [ 

CO(I)=C(I) 

Е (Н(М)) 1,251, 

C(I)2GUND *T C1) *CO( 2) *CO(T) 

GO TO 3 

CONT INUE 

А-ЕХРЕ(-Н(М)%Т(1)) 

8=T (1) *H(N)*HIN)/GIN) 

ClII)=CO(J)%*( (1+2>A)/B-Awx*G(N)/H(N))+C(J)*#*#(G(N)/H(N)-(1s—-A)/B) + 
1CO(I)*A 

CONTINUE 

RETURN 

END ; 

SUBROUTINE LIMITER (NsI»J) 

DIMENSION 6G(99),H(99)5»5P(99)5Q(99) 9R(99) ,C( 99) ,CO(99), COO(99), 
1С000(99)»С0000(99)»Т7(20)»Х1(900)»Х2(900) »Х3(900)»Х84(900)»Х5(900) 
2\1 (900) Ү2(900)%>Ү3(900)»Ү8(900) »Ү5(900) »2(26)»((26)» 
31TITLE(12) s JTITLEC12) SKTITLE(12) s LTITLEC 12) 4 MTITLE(C12)» 
bh X(9550)5»Y(9550) 

СОММОМ б»Н»Р»0»В»С»С0» С00»С000»С0000»7»Х1%9Х2»9Х3»Х4»Х5»Ү1»>Ү2»Ү3» 
1 Ү4»Ү5»2>(»Х>Ү 

С0О000(1)-С000(1) 

СООО(1)=СОО(1) 

С00(1)=С0 (1) 

CO(1)=C(1) 

CLI)2GUIN 2 * C C. 7) 


= AO 
жо 


= шали ра р T E - и ДІ - (ΠΕ Τ᾽ Τη ΠΠ πι p ποτ Γη wm... Texten pt * + tcn 


} 





10 
20 


30 
40 


10 
15 


20 
30 


CHECK UPPER LIMITsTHEN LOWER LIMIT 

IF(C(I)-H(N))10910920 

IF(P(N)-C(1))40+40930 

C(1)=H1N) 

GO TO 40 

C(1)=P(N) 

RETURN 

END : 

SUBROUTINE MAGNITD(N,T»J) 

DIMENSION 6(99)»Н(99)»Р(99)»0(99)К(99) »С(99)»С0(99)»С00(99)» 
1С000(99)С0000(99)»Т(20)»Х1(900)»Х2(900) »Х3(900)»Х4(900)»Х5(900)» 
2Ү1(900)»Ү2(900)»ҮЗ(900)»Ү8(900)Ү5(9001,2(26)»1(26)» 
3ITITLE(12)» CEA o 5 MTITLE(12)» 

¿ X(9,50),Y(9,50) 

СОММОМ СэН»Р» 0, КБ» С»С0»С00»С000%СОООО»Тәхі»Х2»Х3»Х4»Х59Ү1»Ү2»Ү3» 
1 Ү6»Ү5»2>(»Х»Ү 

СО000(1)-СООО(1) 

СООО(1І)-СОО(1) 

COO(I)zCO(I) 

COLI )2C(I) 

C(I)2ABSF(C(J)) 

RETURN 

END 

SUBROUTINE IDRELAY(NoI9J) 

DIMENSION G(99) 9H(99) 9P(99) 9Q(99) 9R(99) C199) »›СО( 99) »СОО (99 } э 
1С000(99) «С0000(99)»Т(20) 9X1( 900) 9X2(900) 9X3(900) »Х4 ( 900} »Х5 (900 } » 
2\1 (900) э\2 (900) э\З (900 ) э\4 (900) 5Y5(900) 9Z(26) si (26) » 
31ITITLE(12) » JTITLEC 120) sKTITLECI2) sLTITLECI2) 4s MTITLE(12)» 

4 Х(9»50)»Ү(9»50) | 

COMMON СОэН»Р» ,0»8»С»С0»С00»С000»С0000»7»Х1»Х2»Х3»Х4»Х5»Ү1»Ү2»ҮЗ» 
1 Ү4»9Ү5,2,1»Х»Ү 

СОООО( 1) «“СОООС(Т) 

СООО( 1) <СОО(1) 

COO(I1)2CO(lI) 

СО(1)=С(1) 

IF(C(J)) 10915920 

SUBSTITUTE LOWER LIMIT 

C( I)zH(N) 

GO TO 30 

C(1)=08 

GO TO 30 

SUBSTITUTE UPPER LIMIT 

C(1)=G(N) ` 

CONTINUE 

RETURN 

END 

SUBROUTINE GAIN(N»I»4J) 

DIMENSION С(99)»›Н( 99} »›Р (99 } ›О (99) эВ (99 } »С( 99 } »›СО( 99) »СОО ( 99 } » 
1С000(99) С0000(99)»Т(20)»Х1(900) »Х2(900) »Х3(900)» Ха(900)»Х5(900), 
2Ү1(900)»Ү2(900) »ҮЗ(900)»Ү6(900)»Ү5(900)»2(26)»((26)» 
31TITLEC12) » JTITLECI12) 3SKTITLECI2) 3LTITLEC12) S MTITLE(C12) » 

4 Х(9»50)»Ү(9»50): 

COMMON GsHsPs3QsRsC»COsCOOsCOOOsCOOOOS T» X1» X295 X25 X64 X5» Y1»Y2» Y23» 

1 Ү4»Ү5,2,1»Х»9Ү 
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CO000( 1)=C0001(1) 

С000(1)-СОО(1) 

COO(1)=CO(T) 

CO(CI)sC(I) 

C(I)=G(N)#C(J) 

RETURN 

END 

SUBROUTINE MULTPLY(N»IsJsK) 

DIMENSION G(99) ,H( 99) »P( 991,QU99) R(99) C( 99), CO( 99) » COO( 99) , 
1С000(99) »С0000(99)»7(201»Х11(900) Х2(900)»Х3( 900)»Х4(900)»Х5(900), 
2Ү1(900) Ү2(900) Ү3(900)»Ү4(900) Ү5(900) »2(26)»0(26)» š 
SITITLE(12) 9 JTITLE( 12) KTITLE(12) sLTITLE(12) oMTITLE( 12)» 

4 X(9+50)9Y(9950) 

COMMON G»HsP ОКС СО СОО COCO E T»X1»X2» X3»X6 »X59»Y1»Y2»Y3» 
1 Ү4,Ү5,2,1»Х»Ү 

СО000(1) «С000(1) 

СО0ОО(1)«СОО(1) 

СОО(1)-СО(1) 

COtI)sCÓI) ` 

C(I)5CCJ) *C(K) 

RETURN 

END | 

SUBROUTINE DIVIDE(Ns» Ir J»K) 

DIMENSION G(99)5H(99) P(99)5Q( 99) 3R(99) C( 99) ,CO( 99) COO(99), 
1С000(99) »С0000(99)»7(20)»Х1(900)»Х2(900) »Х2(900)»Х4(900)»Х5(900)% 
2Y1( 900) 9¥2(900) 9¥3(900)9Y4(900) 9¥5(900) 52(26)5L(26)>» 

BEPRTTLE ( 12)» JTITLECI2 )sKTITLEC(1 2) s LTITLECI2) MTITLE (C120 » 
% Х(9»50)»Ү(9»50): 

COMMON СЕ» С00» СО0ОО» CO000 »T »X1»X2»X3»X4» X5»Y1»Y25»Y23» 
1 Ү8%5Ү5»794»Х 

A 

СООО (1) =С00 (1) | | | 

COO(I)zCOlI) 

CO(I)=C(I) > 

ЕСІ!-СІ.)) /СІК) 

ВЕТУВМ 

ЕМО 

SUBROUTINE CURVE(NOsI9J) 

DIMENSION 6(99)»Н(99)»Р(99)»0(99)»К(99)»С(99)»С0(99) »С00(99)» 
1С000(99) »С0000(99)»7(20)»Х1(900)»Х2(900) »Х3(900)»Х8(900)»Х5(900): 
2Ү1(900) »Ү2 (900) »ҮЗ(900 ) »ҮФ (900 | »Ү5 (900) »2(26) FL (26) 
41TITLEC12) s JTITLEC12) sKTITLEC12) 3 LTITLECI2) sMTITLE(12) > 
4 Х(950)»Ү(9»50)- 

СОММОМ С»Н»Р»СО»ЕВ,»С»С0»С00» СО00»С0000»Т79Х1»9Х2»9Х3;9Х4»9Х5»9Ү19Ү2»Ү3» 
1 Ү4%Ү5»2512»Х>Ү 

С0000(1)-С000(1) 

СО00(1)-С00(1) 

COO(lI)2COLlI) 

CO(1)=C(1) 

THIS-1S.A LINEAR INTERPOLATION USING THE FORMULA 

CU 1)=YO-(YO-YN)(C(J)=X0)/(XN-X0) 

К=0 


15 К=К+1 


IF(51-K) 30, 30,1 
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20 
30 


IF(C(CJ) -X(NO»sK)) 10520515 

IF(K-1) 3053092 

CCI)2Y(NO S K-1)- CYCNOSK-1)-Y (NOS K) ) X (CCJ) -CXCNOS, K-1)) 7 
1 (Х(МОЭК)-Х(МО»К-1)) 

GO TO 6 

PRINT 4»NO 

FORMAT t 30HFUNCTION GENERATOR ERROR CURVE 912) 

PRINT 5»J»CtJ) 4 

FORMAT(3H X(»1353H 12,E20.5) ! 

STOP 

CUI) SYUNO SK ) 

RETURN 

END 

SUBROUTINE RECIP (М,1».) 

DIMENSION G(99) 9H(99) 9P(99) »›О (99) »›В (99 } ›С (99) ›СО (99) ›СОО (99 } » 
1С000 (99) ›СОООО0 (99) ›»Т{20)»›Х] (900) ›Х2 (900) ›Х 3 (900) ›Х4 (900) ›Х5 (900 } › 
2Ү1(900) »Ү2(900) »ҮЗ (900 ) »Ү4 (900) »Ү5 (900) »2 (26) »_ (26) » 

3 ITITLE(12)» JTITLEC12) »KTITLE(12) ;LTITLE(CI2)  MTITLE(12), 
4& X(99550)5Y(95,50] 

COMMON Gs»H5»P»Q$* R» C*CO»COO»s»COOO»COOOOS T»X1»X2»X23»X^49 X5» Y1»Y2» Y3» 
1 Ү4»Ү5,79(»Х»эҮ 

COOOO0(I)zCOOOlI) 

COOO(I)sCOO(lI) 

COO(lI)zCO(lI) 

СО(1)=С(1) | 

C(I)=G(N)7/C (4) | 

RETURN Р 

ЕМО 

SUBROUTINE РЕ ΑΥ(ΝΟ9Ι 9) 

DIMENSION 6(99)»Н(99)»Р(99)»0(99),ЕК(99)»С(99)»С0(99)»С00(99)» 
1С000(99} »›СООО0 ( 99) ›Т( 20} + Х1 (900) »Х2 (900) 9 Х3 (900) »Х4 (900) »X5( 900) > 
2\1 (900) »\2( 900 ) »\З {900 ) э\4 (900 } 9¥5( 900) 92(26) oh (26) 9 
31ITITLE(12)»JTITLE(12) »KTITLE(12)»LTITLE(C12) 3MTITLE( 12)» 

4 Х(9»950)»Ү(9»>50) 

СОММОМ б»Н»Р»0»В»С»С0>С00»С000»С0000>79Х1»Х2»Х3»Х4»Х5»9Ү1»Ү2» УЗ» 
1 Ү6»Ү5,)»7,9(»Х»Ү 

С0000(1)-С000(1) 

COOO(I)=COO (I) 

COO(I)sCOlI) 5 

СО(1)=С(1) 

IF(C(99)-GIN)) 10920920 

C{1)=0.0 

GO TO 30 

C(I)=C(J) 

CONTINUE 

RETURN 

END 

SUBROUTINE ADDER(N»I»JsK»LK) 

DIMENSION G(99)»H1(99)»P(99)90Q199)»R199)9C199)»C0199) »C00199)» 
1С000(99)»С0000(99)»Т (20) »Х1(900)»Х2(900) »Х3(900)»Х4(900)»Х5(900)» 
2Ү1(900)>Ү2(900)»Ү2(900)»Ү6(900)»Ү5(900) »2(26)>1(26)» 
31TITLE(12)» JTITLEC 12) KTETLE (12) oLTITLE(12) sMTITLE( 12)» 


4 Х(9,50)»Ү(9»50): 
СОММОМ б»Нн»Р»0»В»С>С0»С00»С000»С0000»Т»Х1»Х2»Х3;%Х4» Х5%9 Ү1»Ү2»Ү3» 
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1 Ү4,Ү5,2,(»Х»9Ү 

СОООО(І)-С000(1) 

СООО (Т) =СО00 (1) 

COO(I)=CO(I) 

CO(I)=C(I1) 

IF(J) 19393 

JJ=XABSF( J) 

S1=-C (JJ) 

GO TO 4 

S1=C(J)' 

IF(K) 59696 

KK=XABSF(K) 

S2=~-C( KK) 

GO TO 7 

52-С(К) 

IF(LK)859,10 

LLzXABSF(LK) 

$3z-C(LL) 
GO TO 11 

S3zC(LK) 

C(I)=S1+S2+S3 

GO TO 15 

C(1)=S1+S2 

CONTINUE 

RETURN 

END 

SUBROUTINE DELAY1L(NoI9J) 

DIMENSION G(99)5,H(99)5P(99)5Q( 99) 9R( 99), C( 99),CO( 99) COO(99) » 
1С000(99)»С0000(99)»Т(20)»Х1(900) »Х2(900)»Х3(900)» Х4(900)»Х5(900)» 
2Ү1(900)»Ү2(900)»ҮЗ3(900)»Ү4(900)»Ү5(900)»2(26)>((26)» 
З1ТТТЕЕ (12) эУТТТЕЕ (12) эЭКТ1ТЬЕ (12) 5LTITLE(12) 4$ MTITLE(12)» 
à X(9550)»Y(99550) »Al3»1002) »M(5} | 

COMMON G»HsPsQsR«»C»COSsSCOOsCOOO»sCOOQOOS T» X1» X2»X35» X69» X5» Y1»Y29»Y3» 
1 Ү4»Ү5,7,(»Х»Ү 

ІҒ (С(99)-Т(1)) 291,1 

DO 5 KK=193 

DO 5 JJ=191001 

A(KKeJJ)=C(1) 

НМ) =6(№) /7(1) +20 , 
MIN) H(N) 

CO000( 1)=C000 (1) 

CO0O0011)=C00 (1) 

СОО(1)=СО(1) 

СО(1)=С(1) 

MM=M(N) 
AUN»MM)=C(J) 
DO 3 K=1»MM 

AUNS K ) SÀUN S K £1) 

C(I)=A(N>1) 

, RETURN 

END 

SUBROUTINE NOISE(N»I) 

DIMENSION G(99)5H(99) »P(99)50(99) 9R(99),C(99),CO(99) ,COO( 99), 
1С000(99)»С0000(99) Т(20)»Х1(900) »Х2(900)»Х3(900)»Х4(900)»Х5(900)» 





no 
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2Y1(900) »Y2(900) » Y2 
3ITITLEC12) »JTITLE( 
4& X(95,50)5»Y(95,50) 

COMMON 65»Н»Р»0»К»С»С0»С00»С000»С0000»7»Х1»Х2»Х3»Х8»9Х5»9Ү1»>Ү2»Ү3» 
1 Ү6»Ү5»291»Х»Ү 

С0000(1)-С000(1) 

COOO(I)2COO(I) 

COO(I)2CO(I) 

CO(I)=C(1I) 

CALL RNDEV( 1220703125 sDEV) 

C(I) =H(N) *DEV+GI(N) 

RETURN 

END 

SUBROUTINE ERROR (NoKosIo0J) 

DIMENSION G(99),H(99),P(99),Q(99) 9R(99),C(99),CO(99),COO(99), 
1С000(99)»С0000 (99) ›Т(20) ›Х1 (900 } »Х2 (900) »ХЗ (900) ›Х4 (900 } »›Х5 (900 )} » 
2ү1 (900) »Ү2(900) »ҮЗ (900 ) »Ү4 (900) »Ү5 (900) »2 (26) » (26) » 
ZITITLE(12) sJTITLE(12) sKTITLE(12) sLTITLE(12) sMTITLE(12) 5 
4 Х(9;50)»Ү(9»50) 

СОММОМ С»Н»Р»0»В»С»С0»С00»С000»С0000»7»9Х1»Х2»Х3»9Х4»Х5»Ү1>Ү2»Ү3» 
1 Ү4»Ү5,7,(»Х»9Ү 

THIS ROUTINE SUMS(C(I)-C(J))**2(DT) OVER T AND DIVIDES BY T 

IF(C(99)-T(1)) 25152 

ERR=0e0 

ERRSERR*(UCCI)-CUU) ) **2) *T (1) 

C(N)=ERR/C(99) 
C(K)=C(1)=C(J) 

IF(C(99)-T(2))59393 

FORMAT (1H1) 

PRINT 6 

ERX=ERR/T (2) 

FORMAT(11HTHE NUMBER 912910H ERROR IS »Ғ20.9) 

PRINT 4 sNsERX 

RETURN 

END 

SUBROUTINE CMPLXPL (NoIoJ) 

DIMENSION G(99),H(99),P(99)5Q( 99) 9R(99) Cl 99) , CO( 99) COO( 99) , 
1С000(99)»С0000(99)»Т7(20)»>Х1(900) »Х2(900) Х3(900)»Х4(900)»Х5(900)» 
2Ү1(900)»Ү2(900)»Ү2(900)»>Ү4(900)»Ү5(900)»2(26)»>((26)» 
31TITLE(12)>JTITLE(12)»KTITLE(12)>»LTITLE(12)»MTITLE(12)» 

4 X(9+50)»Y(9950)»AK(492)9YC(2)9DY(2) 

СОММОМ б»Н»Р»0»В»С»С0»С00» С000» С0О00О»Т» Х1>Х2»Х3»Х6%Х5»Ү1>Ү2>Ү3» 
1 Ү8»Ү5»29(»Х»Ү р 

СО000 (1)-С000(1) 

COOOlI)2sCOO(I) 

СОО(1) «СО(1) 

CO(I)=C(I) 

Y5(700*2*N) x6UN)*CO( J) -HUN) *X5( 700*2X*XN) -P(CN) *CO(T) 

X5(7004-2*N)2Y5( 700-2*N)*TC1) *X5( 7004 2*N) 

C(1)=X5(700+2*N)*T(1)+C0(1) 

RETURN 

END 
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Appendix II 


Root Locus Program 
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PROGRAM RLOCUS 
PROGRAMMERS RM NUTTING » JL FENICK » AND JA POPE 
THIS PROGRAM WILL PLOT A ROOT LOCUS FOR A CHARACTERISTIC EQUATION UP TO 
ORDER 30. ROQT LOCUS POLES ARE PLOTTED WITH AN X» ROOT LOCUS ZEROS ARE 
PLOTTED WITH A SQUARE» AND INTERMEDIATE ROOT POINTS ARE PLOTTED WITH A 
PLUSe THE STARTING VALUE OF ROOT LOCUS GAIN AND THE NUMBER OF DECADES 
TO BE SPANNED BY THE GAIN MUST BE SPECIFIED. THE GRAPH PLOT IS BASED ON 
PLOTTING EVERY TENTH POINT AS THE GAIN VARIES BETWEEN ITS INITIAL AND 
FINAL VALUE IN 300 STEPS. 

THE DATA CARDS ARE SUBMITTED IN THE FOLLOWING MANNERe SUBMIT A 
COMPLETE SET OF DATA CARDS FOR EACH ROOT LOCUS TO BE PLOTTED. 
CARD 1 THE FIRST LINE OF THE GRAPH TITLE (IN COLUMNS 1-48) 
CARD 2 THE SECOND LINE OF THE GRAPH TITLE (IN COLUMNS 1-48) 
CARD 3 THE ORDER OF THE CHARACTERISTIC EQUATION (12 FORMAT) 
CARD 4 CONSTANT COEFFICIENTS IN DESCENDING ORDER. (8Е10.5 FORMAT) 
CARD 5 COEFFICIENTS OF THE VARIABLE IN DESCENDING ORDER (8E10,5 FORMAT) 
CARD 6 INITIAL VALUE OF THE VARIABLE ( E10,5 FORMAT) » MUST NOT BE ZERO 
CARD 7 NUMBER OF DECADES TO BE SPANNEDe (FROM 1-10) » (12 FORMAT) 
CARD 8 GRAPH SCALE TO ONE SIGNIFICANT FIGURE». (£1065 FORMAT) 
DIMENSION R(129)5X(129)» IT(10) *ROOTR(128) sROOTI (128) » ITITLE( 12)» 


1A(129)5B( 129) sROOTJ( 128) sROOTM( 128] 


15 
206 


14 
200 


204 


24 


- 
> 


28 


205 
207 
22 


23 


COMMON R»VAR»NO+ROOTR»ROOT 1 
DO 15 K=19129 
MOD=1 
LAB=4H 
MM=0 | 

DO 14 L=1»10 

IT=8H 

FORMAT (6A8 ) 

FORMAT (12) 

FORMAT (E10-5) 

ВЕАО 2009 (ІТІТІЕ(І)»іІг1»6) 
ВӘЕАр 2009 (ІТІТІ-Е(І)» Із7912) ; 

FORMAT (1H19»////»317HTHE INPUT DATA ISs»////) 

PRINT 24 

FORMAT (///»36HORDER OF ΤΗΕ CHARACTERISTIC EQUATION»///) 
PRINT 28 

READ 203»МО | 

PRINT 203,NO | 

‚ М=МО+1 

FORMAT (8E10. 5) 

FORMAT (8E12.5) 

FORMAT(///s41HCONSTANT COEFFICIENTS IN DESCENDING ORDER»///) 
PRINT 22 

READ 205.5 (A(K) »Kz 1o N) 


PRINT 207»(A(K) »Kz15sN) 
FORMAT (/// » &BHCOEFFICIENTS OF THE VARIABLE IN DESCENDING ORDER» 
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25 


26 


27 


201 


202 


10 


ή] 


PRINT 23 

READ 2059(B(K) sK=19N} 

PRINT 2079(B(K) sK=19N) 

FORMAT(///»29HINITIAL VALUE OF ΤΗΕ VARIABLE>///) 
PRINT 25 

READ 204>VAR 

PRINT 204,VAR 

FORMAT(///»31HNUMBER OF DECADES TO BE 5РАММЕ0»///) 
PRINT 26: 

READ 203,ND 

PRINT 203»ND 

ЕОНМАТ(///»5Н5СА(Е»///) 

PRINT 27 ` 

READ 204,5 XSCALE 

PRINT 204,XSCALE 

YSCALE=XSCALE 

FORMAT(21HITHE SYSTEM POLES ARE ////) 

PRINT 201 

CALL ROOTS2 (AsXsNOs» ITsROOTRSROOTIS,MM»s—-.55*.5) 


CALL DRAW(NO,ROOTR,ROOTI»MOD,1»LAB, ITITLE»,XSCALE,YSCALE» 


11»6»2»2»17»8,1»1 А5Т) 


MOD=2 

FORMAT(//»21H THE SYSTEM ZEROS ARE»///) 
Қ-1 

IF(B(K)) 1921 

K=K+1 7 

GO TO 3 

NORD=N-K 

IF(NORD-1) 65455 

ZERO=-B(K+1)/B(K) 

FORMAT (//»16HTHE SYSTEM ZERO=9+E100+5»///) 
PRINT 7»ZERO 

GO TO 8 

PRINT 9 

FORMAT(//s25HALL ZEROS ARE AT INFINITY) 
GO TO 8 

NN=NORD+1 

DO 10 L=1»NN 

R(L)=B(K) 

K=K+1 

PRINT 202 

СА ВООТ52 (В»Х»МОНКОр» ІТ»КООТМ»КО00Т7.»ММ»-45945) 


CALL DRAW{INORD»ROOTM,»ROOTJ»MOD,3,LAB,ITITLE,XSCALE,YSCALE» 


11»6»2»2»7s8»1»LAST) 


|. MOD=2 


8 
31 
32 


CONT INUE 
GO TO(31»32533,345,35536»375,38539,40) ND 
G=1.0076 
GO TO 41 
6=1.016 
GO Т0 41 


33 6=]1.0245 


60 ТО 41 


——————— = -- — ἜΤ г М: DOING тг 
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34 
35 
36 
37 
38 
39 
40 


41 
30 


42 


60 


300 
50 


69 


70 
71 


100 
29 


98 
101 


20 


G=1+0312 

GO TO 41 

G=1.0394 

60 ТО 41 

6=1.0483 

GO TO 41 

621640568 

GO TO 41 i 
G=1.0633 ` 
GO TO 41 

б=1,6071 

GO ТО 41 ` 

G=1.078 

PRINT 30 

FORMAT (1H19////s61HROOTS FOR THE SPECIFIED VALUES OF THE VARIABLE 


LARE AS FOLLOWS, 1117) 


PRINT 42 | 
ΕΟΒΜΑΤ(1ΟΧ»ώΗ УАК»4Х»ЭНКЕА( РАВТ»3Х»9НІМАС РАЕВТ» 


113X»9HREAL PART»>3X»9HIMAG PART>13X»9HREAL PART»>3X»>9HIMAG РАКТ»///) 


DO 101 J=1»30 

DO 100 K=1510 

FORMAT(1PE12¢3) 

PRINT 60,VAR 

DO 300 L=1»N 

R(L)=A(L)+B(L)VAR 
FORMAT(15X>3(1P2E12.3»10X)) 

CALL ROOTX 

DO 71 JJz1»NO 
IF(ABSF(ROOTI(JJ))-54.,E-0^) 69.69.70 
ROOTJ(JJ)s0, 

GO TO 71 

ROOTJ(JJ)sROOTI(JJ) 

CONTINUE 

PRINT 50»(ROOTR(I)»ROOTJ(1)>1=19NO) 
VAR=G*VAR 

IF(J-1) 1015101599 

MOD=2 | 
IF(J-30) 101598598 

MOD=3 

CALL DRAWINO»ROOTR»ROOTJ»MOD»2»LAB» ITITLE+XSCALE» YSCALE» 


1: 1»6»2»2»7»8»1»LAST) 


60 ТО 206 

ЕМО 

SUBROUTINE ROOTX 

DIMENSION C(31)50(29)5R(129) *ROOTR(128) »ROOTI (128) 
COMMON R»VAR»NO»ROOTR»ROOTI 
M=1 

BETAN sROOTI(M) %.001 
ALFAN=ROOTR(M) 

DO 7 1=1»100 

Sz2.*ALFAN 
T=-tALFAN##*2+BETAN*#2) 
C(1)=R(1) 

C(2)=R(2)+S*R(1) 


1 
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17 


21 


18 


19 


οι 
Eg o чм + 


13 


12 
15 


16 


NC=NO+1 
DO 2 L=3,NC 
C(L) = R(L)+S*C(L-— 


BN= BETANXCINO) 
IF (NO-3) 21» 17 


1) 


18 


| * T*C(L-2) 
AN= C(NO+1)-ALFAN*CINO) 


СМ = 3,.*RUL)*(ALFAN**2-BETAN**2) + 2,.*R(2)*ALFAN * R(3) 


ОМ = 6,.*R(1)*ALFAN*BETAN * 24,*R(2)*BETAN 


DO 3 N=3,NU 


D(N)=C(N)+S*D(N-1)+T*D(N—2 ) 


GO TO 19 

CN = 2e*R(1)*ALFAN + R(2) 
DN = 2e*#R(1)*BETAN 

GO TO 19 

0(1) - С(1) 

D(2)2C(2) *S*D(1) 

NU=NO- 1 


CN= С(М№О)-2.*0( м0-2 ) *ВЕТАМ* *2 


DN= 2.*BETAN*(D(NO-1)-ALFAN*D(NO-2)) 
ALFAsALFAN-(AN*CN«-BN*DN)/(CN**240Nx**2 | 
BE TA=BETAN+( AN®DN-BN#CN ) / (CN¥#2+DN##2 ) 
IF (ABSF(ALFA-ALFAN)-5.E-4) 49495 
IF(ABSF(BETA-BETAN)-5.E-4) 


ALFAN=ALFA 
BETAN=BETA 
PRINT 50 


FORMAT (46H NO CONV 


GO ТО 12 
ROOTR(M)=ALFA 
ROOTI(M)=BETA 


IF (ABSF(ROOTI(M))-5.E-4) 12912913 
ROOTR(M+1) = ROOTR(M) 
ROOTI(M+1) = -ROOTI(M) 


М= М+ 1 


IF(M-NO) 15516516 , 


М=М+ 1 


| 


|| 


` 


C 


GO TO 20 
RETURN 
END 

END 
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Appendix III 


Subroutines INVAL and DIAGRAM and Data for Cruising 
Conditions and Ninety Per Cent Delta Plant 
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SUBROUTINE INVAL 


ОІМЕМОІОМ 6(99)»Н(99)»Р(99)»0(99)В(99)»С(99)»С0(99)»С00(99)» 
1С000(99)»С0000(99) ЭТ(10)»Х1(900)»Х2(900) »Х3(900), 
2Ү1(900)»Ү2(900)»Ү3(900)» 
3ITITLE(12)sJTITLE(12)yKTITLE(12)) 
4 Х(950)»Ү(9»50) 

COMMON, бЭН»Р»0»ЕВ»С»СО»С00»С000»С0000»79Х1»Х2%Х3%9Х»9Ү»Ү1»Ү2%ҮЗ 

С(50) ES READ FROM IC CARD 

C(42)=159006 
C(97)215. 

С(92)-С(50) 

CALL CURVE (2»47,50) 

CALL CURVE (1»76,»50) 

CALL ADDER (48594859 769239-47) 

C(77)=0.04*C(48)-3360 

C(49)=C 177) 

C(52)=C197)-C1(49) 

С(531-(6(53)/Н(523))%С(52) 

С(78)-0.0001142 %С(50) 43.0 

C(57)2C(78) 

C(54)2C(53)4C(57) 

С(55)-С(54) 

C(67)=C(55) 

C(66)=C(67)-3.0 

С(3)-С(66)%(Н(66)/6(66)) 

C(65)=C(3) M 

С(79)-С(65) 

С(75)-С(79)/.254 

C(72)2C(75)/G( 75) 

C(56)zC(55)-C(67). 

С(703-С(72)%(Н(72)/6(72)) 

C(69)2C(70) 

С(63)-С(69)%(Н(69)/6(69)) 

C(61)=C(63) 

C(60)=C(61) 

C(58)2C(56)*(G( 58) /P( 58)) 

C(51)2C(58)*G(51) 

С(59)- С(60)-С(51) 

ІҒ(С(55)-С(67))191»2 , 

С( 64) =С (55) 

GO TO 3 

C(64)=C(67) 

CONTINUE 

FUEL OIL FLOW (GF) FOR CRUISING SPEED «4140. 

C(46)=C(64)-300 

C(45)=C (46) *H(46) /G(46) 

c(39)=C(45)/0.00141 

С(81)-С(42)-6(40)%С(39) 

C(40)2C(39) *G(^0) 

C(90)2C(239) 

С(38)-С(41)-С(39) 

C(87)2C(38) 

C(80)= C(38)*0,00141 +30 

C(81)=0200141%*C(41)4+360 

C(44)=C (81) 

C(43)=C(80) 

C(5)=C(14)+1500 

C(33)=C(5) * 

C(93)=C150)-" 

C(4)=0.000114*C(93)+300 





С({29)=С(4) 
С(30)-С(29) 
С(31)-С(57)-С(30) 
C(32)=C(31)*(G(32)/H(32)) 
ае Са 01:93) 
C(25)=C(24)*(G(25)/P(25)) 
RETURN 

END 
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1000 
2000 


3000 
4000 
1100 
2100 
3100 
4100 
1200 
2200 


3200 
4200 


CALL REALPL (6656653) 


| 


SUBROUTINE DIAGRAM 
DIMENSION 6(99) эН(99) Р(99)»0(99)»К(99)»С(99)»С0(99)»С00(99)» 
1С000(99) » С0000(99) »Т(10)»9Х1(900)»Х2(900)»Х3(900)» 
2Ү1(900) «Ү2(900)»ҮЗ(900)» 
ZITITLE( 42) sJTITLE(12) sKTITLE(12)>5 
4 Х(9,50)»Ү(9»50) 
СОММОМ бСэН»Р» 08» С»С0»С00» С0О00»С0000»Т7»9Х1»Х2»Х3%9Х»9 Ү» Ү1»Ү2»Ү3 
ІҒ(С(99)-5.0)10910»20 | 
С(50)=56000.- (600. *С (99) ) 
CONTINUE | 
CALL CURVE (2947,50) 
CALL CURVE (1976550) : 
CALL ADDER (485485,765»235-—^7) 
СО(77)-С(77) 
С(77)-.04%ҰС(48)-33,0 
CALL CMPLXPL (49 949577) 
С(97)-15. 
CALL ADDER (525525,975,-49,0) 
CALL REALPL (5295329525) 
CO(78)2C(78) 
С(78)-.0001142%С(50)%3.0 
CALL CMPLXPL(57957978) 
CALL ADDER (5455455355730) 
CALL CMPLXPL (55555554) 
CALL ADDER (563569555-6790) 
CALL CMPLXPL (58,58,56) 
CALL GAIN (51551558) 
IF(C(51)2-30.,0 )1000»91000»2000 
IF(C(51)-0. гооо» док» 4000 
C(51)=30.0 | 
GO TO 4000 E 
C(51)20,0 
CONTINUE ч 
CALL REALPL (59959958) 
IF(C(59)-30,0)11005»110052100 
ІЕ( С(59)-0.0)3100»3100»4100 
С(59) = 30,0 , 
GO TO 4100 
С(59)-0.0 
CONTINUE 
CALL ADDER(60560,5155950) 
IF (C(60)-30.) 1200»1200»2200 
IFtC(60)-0.) 3200»3200»4200 
С(60)-30. | 
СО ТО 4200 
C(60)=06 
CONTINUE 
CALL CMPLXPL (61961560) 
CALL DERIVA (62962961) 
CALL ADDER (6356356156250) 
THE AIR SYSTEM GOES HERE 
CALL REALPL (69369963) 
CALL CMPLXPL (70570569) 
CALL REALPL (72972970) 
CALL GAIN (75975972) 
СО (79) =С (79) 
C(79)=C(75) %.254 
CALL REALPL(65565579) | 
CALL DELAY1 (3»3»65) | 


| 





C(67)=C(66) +3.0 
IF(C(55)-C(67))1»1»2 
1 C(64)=C(55) 
GO TO 3 
2 C(64)=C(67) 
3 CONTINUE 
С OIL SYSTEM GOES HERE 
CO(34)=C(34) 
C(34)zCÍ64)-(C( 46) 43.0) 
CALL REALPL (35 935534) 
CALL GAIN (36936534) 
CALL ADDER (3793793593690) 
CALL REALPL (865986537) 
CALL ADDER( 38389-86987 ,50) 
CALL ADDER( 39 939 9-38 94150) 
CALL ADDER(89,89,39,-90,5,0) 
CALL GAIN(409+40939) 
CALL ADDER( 41941 5-4094250) 
СО(81)-С(81) 
СО (80) =С (80) 
C(80)= C(38)*0.00141 +3.0 
С(81)-0.00141%С(41)-3.0 
CALL CMPLXPL(43 943580) 
CALL CMPLXPL(44,44,8]) 
CALL ADDER (459+459+44»>-4390) 
CALL REALPL (46546545) 
C WATER SYSTEM 
CO(5)zC(5) 
C(5)=C(14)+15.0 
CALL CMPLXPL (3393395) 
CALL ADDER (249249=-33932997) 
CALL CMPLXPL( 84984 524) 
CALL GAIN(25925984) 
CALL REALPL( 26926984) | 
CALL ADDER (2792792592630) 
CALL DELAY1 (292927) 
CALL CMPLXPL(2892892) 
CO(4)=C (4) 
С(4)-0.000114%(С(93) «С(28))%3.0 
CALL CMPLXPL(29 929 94) 
CALL REALPL (30930929) 
CALL ADDER(315»315»575»-3050) 
CALL REALPL (2292293291) 
С WATER LEVEL 
CALL ADDER(915915505,-9250) 
CALL REALPL(696991) 
CALL REALPL( 797991) 
CALL CMPLXPL(85857) 
CALL ADDER (99998s-690) 
CALL REALPL (10510528) 
CALL CMPLXPL(11»115»10) 
CALL REALPL( 12912989) 
CALL CMPLXPL(13913912) 
CALL DELAY1 (151513) 
CALL ADDER (145»1459511»1) 
C STEAM PRESS 
CALL REALPL (83983991) 
CALL REALPL (15.15.83) 
CALL REALPL (16916915) 
CALL ADDER (17»917;15»-16>0) 
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CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


REALPL (82582589) 

REALPL (18918582) 

REALPL (19519518) 

ADDER (20»20»18»-19»0) 
REALPL (21521528) 

REALPL (22922921) 

ADDER (23,23,20,-17.-22) 


RETURN 


END 





E 


DATA FOR CRUISING CONDITIONS DELTA PLANT 


-48 

01 0.8 

02 025 

Ш 0.5045 

06 1.12 ES 

2M .3805 -05 0.3141 

ОЧ 1.0 0.763 

10 5.54- 07 

ЕМ 1.0 0.422 

12 1.86 -04 2.22 

13 1.0 0.76 

15 3.68 -04 .0445 

ШОР 1.2282 1.41 

18 7.87 -03 0445 

19 1.732 1.76 

21 ie 13 -06 

22 1.0 0.0445 

571.26 

26 0.1 

28 9360. „531 

19.2 26.2 

30 0.22 0.222 

32 14„7 29. 4 

33 4.0 204 

2570.392 

36 2.82 

40 0.0445 

43 30.3 5.9 

44 30.3 5.9 

46 30.87 29.4 

49 50.0 «6 

51 .944 

53 73.6 29.4 

01250 e 9.5 

57 19.25 26.15 

58 333,3 1363 

59 «106 

61 100, 7.9 

62 4.89 

65 25 2.5 

66 30.87 29.4 

69 512.0 1.0 

70 2.0 0.56 

72 000683 el 

75 16134 

82 1.0 

83 1.0 

84 1.25 0.675 

86 24804. 31.2 

2 МОМВЕК ОР 

09 ВКЕАК PTS 
0.0 m 
34000. 


56000. 


NUMBER OF BLOCK DATA CARDS 


1.588 
"0.495 


0.4 


e 7/81 
19, 2 


30.3 
303 


50. 
250. 
19, 25 

D9363 


100. 


1.25 


CURVES 


BOILER PRESSURE 


1200, 
1203, 
1208, 
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CURVE 1 





PRESSURE DROP 


66500, 11; 
100000. 1226. 
120000, 12275 
140000. 1250. 
152000. 1259, 
167000, ολ... 
44 BREAKPOINTS SUPERHEATER 
0.0 DO 
20000. ө О 
28000. 0 
34000. ο 0 
39000. 4.0 
45000. De 
48000. 6.0 
50000. 6.5 
52000. 7.0 
58000. ты 
56000. 8.0 
58000, 8.6 
60000.» 9.2 
62000. 10.09 
64000. 10.5 
66000. IA 
68000. ПИ 70 
70000, 12, 8 
74000. 15420 
76000. 15.0 
78000» 15.9 
80000. 16.6 
84000. 18.0 
86000, 19,0 
88000, 20.0 
90000 ο 21.0 
94000. 23.0 
960006 28.0 
98000. 25.0 
100000, 26.0 
108000. 30.0 
117000. 950 
120000.» Bor ui 
125000. 80.0 
130000, 44.0 
140000. 50.0 
145000. 54.0 
150000. 59.0 
160000. 66.0 
165000. 70.0 
170000. 74.0 
180000. 84.0 
190000. 92.0 
200000» 101.0 
-01 
50 56000. 
0.01 ТІМЕ ІМСКЕМЕМТ 
55.0 TOTAL TIME 
le 
0.2 PLOT TIME 
0193 VARIABLES TO BE PRINTED 


03 OBRIEN GH DELTA SYSTEM LINEAR AIR 
AIR FLOW (C-75) (PERCENT) VS TIME 


9975 


P 


141 


CURVE 2 


(МЕС КАМР ) 3-17 





29897 DELTA F - H 
UEL FLO 
9928 и С-В TALEB 
DELTA FEEDWATER FLOW (c-28) (LBS/HR) VS TIME? 
IME 


142 





10 


1000 
2000 


3000 
4000 
1100 
2100 
3100 
4100 
1200 
2200 


3200 
4200 


SUBROUTINE DIAGRAM 
DIMENSION G(99)5»H(99)»P(99)5Q( 99) R( 99) C( 99)5»CO( 99) » COO( 99) » 


1С000(99) ›СООО0 (99) ›Т(10) ›Х1 (900 } »›Х2 (900) ›ХЗ( 900) » 
2\1 (900) »у2 (900) ›УЗ (900 } »› 


MITITLEELI2YSJTITLE(I2)»sK TI TEE 2 ss 
46 Х(9»50)»Ү(9»50) 
COMMON Gs»Hs»P5»Q5»R»C»CO5»COO»COOO»COOOO» T» X1» X2» X29» X»5Y»Y1»Y2»Y3 
IF(Ct99)-5.0)10,10»20 
С(50)ч152000. -(1600. % С(99)) 
CONTINUE 
CALL CURVE (2947550) 
CALL CURVE (1576550) 
CALL ADDER (48,48, 76»23,-47) 
СО(77)-С(77) 
С(77)-.04%С(48)-33.0 
CALL CMPLXPL(49,49,77) 
C(97)=158 
CALL ADDER (529529979-4930) 
CALL REALPL (53553552) 
СО (78) =С (78) 
С(78).0001142%С(50)-23.0 
CALL CMPLXPL (57957578) 
CALL ADDER (5455455355750) 
CALL CMPLXPL (55555554) 
CALL ADDER (5655695596790) 
CALL CMPLXPL (585958556) 
GALL GAIN (51551558) 
с (51)-30.0 )10005100052000 
IF(C(51)-040)30005,300054000 
8051) =30.0 
GO TO 4000 
€(51)=0.0 
CONTINUE 
CALL REALPL (59959958) 
ІЕ(С(59)-30.0)1100»91100»2100 
ІЕ(С(59)-0.0)3100»3100»4100 
С(59)-20.0 
50 Т0 4100 | : 
С(59)-0.0 
CONTINUE 
CALL ADDER(60»6055155950) 
IF (C(60)-30.) 12005120052200 
ІҒ(С(60)-0.) 2200%23200»8200 
С(60)-30. 
GO ТО 4200 
С(60)-0. 
CONTINUE 
CALL CMPLXPL (61561560) 
CALL DERIVA (62562561) 
CALL ADDER (6396396196290) 
THE AIR SYSTEM GOES HERE 
CALL REALPL (69969963) 
CALL CMPLXPL (70970969) 
CALL REALPL (72572570) 
CALL GAIN (75975972) 
C(9&)=C(75)%1 2 
CO(79)=C(T9') 
C(79)2C( 75) *4,254^ 
CALL REALPL(65»65» 79) 
CALL DELAY1 (3523565) 





CALL REALPL (66566523) 
C(67)=C(66) 43.0 
IFIC(55)=-C(6 71) l7 152 
1 CV64) 2655) 
GO TO 3 
2 C(6&)=C(67) 
3 CONTINUE 
(С OIL SYSTEM GOES HERE 
CO(34)=C(34) 
C(34)=C(64)-(C(46)4+3.0) 
CALL REALPL (35935534) 
CALL GAIN (36936534) 
CALL ADDER (3793733553630) 
CALL REALPL (865865937) 
CALL ADDER(38538,—-865,87,0) 
CALL ADDER( 39,39,—-38,4150) 
CALL ADDER(89589,39,—-90,0) 
CALL GAIN(40 540539) 
CALL ADDER( 41 941 9-4094250) 
CO(81)=C(81) 
СО (80) =С (80) 
C(80)= C(38)*0.00141 -3.0 
С(81)-0.00141%С(41)%3.0 
CALL CMPLXPL(43,43»80) 
CALL CMPLXPL(44 544981) 
CALL ADDER (45945544 5~4350) 
CALL REALPL (46946545) 
C WATER SYSTEM 
СО(5)=С(5) 
C(5)=C(14)+15.0 
CALL CMPLXPL (33593395) 
CALL ADDER (24»24»-33»32»97) 
CALL CMPLXPL(84584924) 
CALL GAIN(25525 584) 
CALL REALPL (26526584) 
CALL ADDER (2792792592630) 
CALL DELAY1 (2592527) 
CALL CMPLXPL(2852852) 
CO(4)=C(4) 
C(4)=0.000114*(C(93)+C(28))+3.0 
CALL CMPLXPL(29 92994) 
CALL REALPL (30530529) 
CALL ADDER(319319579-3090) 
CALL REALPL (32932931) 
С WATER LEVEL 
CALL ADDER(919915509-92 50) 
CALL REALPL (656591) 
CALL REALPL (7537591) 
CALL CMPLXPL(85897) 
CALL ADDER (998-690) 
CALL REALPL (10510528) 
CALL CMPLXPL(115»115»10) 
CALL REALPL (12512589) 
CALL CMPLXPL(135123512) 
CALL ОЕГАҮ1 (1›,1›,13) 
CALL АОБОЕК (14,14›,9,11›,1) 
€ STEAM PRESS | 
CALL REALPL(83583,91) 
CALL CMPLXPL (15.15.91) 
CALL REALPL (16516515) 





CALL ADDER (1751751591690) 
CALL REALPL(82 582589) 

CALL REALPL (18518582) 
CALL REALPL (19319918) 

CALL ADDER (20320918 3-190) 
CALL REALPL (21521528) 

CALL CMPLXPL (22922321) 
CALL ADDER (23923920s-179-22) 
RETURN ` 
END 
END 
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DATA FOR 90 PER CENT FULL POWER DELTA PLANT 


-48 NUMBER OF BLOCK DATA CARDS 

01 e 26 

02 0.25 

03 0,045 

06 1.26 E-06 

07 1.726E-04 0.55 

08 1.0 1.645 1.69 

10 2.0 Е-06 

ШІ 1.0 «508 159 

ШОО 3 {5Е-О05 0.5 

13 10 0.5 „25 

15 .00005 .882 2396 

16 1.985 

18 .00806 .4 

19 1.188 1.41 

21 1.0 E=07 

e2 1.0 0.06 0.01 

#5 1.26 

26 0,1 

28 9360. „531 ο 81 

29 19.2 2.62 19% 2 

2020.222 0.222 

32 14.7 29.4 

33 4.0 2.4 4.0 

350.392 

36 2.82 

40 0.0445 

43 30.3 9 φ 9 30.3 

44 30,3 De 30.3 

46 30.87 29.4 

49 50.0 „б 50. 

51 «944 

М2 13.6 29.4 

55 2506 9.5 250. 

57 1902 2.62 19.2 

255 33363 13.3 25334145 

59 „106 : 

61 100. 7.2.9 100. 

62 „89 

65 25 25 

66 30.87 29.4 

69 950. 1.0 

70 2.0 0.56 "PSU 

72 2.000865 0.1575 

75 1.134 

82 1.0 

83 IRO 

84 1.25 0.675 Ls 22 

86 24804. 31.2 

2 NUMBER OF CURVES 

09 BREAK. PTS BOILER PRESSURE 
0.0 a 1200» 
34000. 1203. 
56000. 1208. 


CURVE 1 





PRESSURE DROP 


66500. ιτ... 
100000. 1226. 
140000, \ 12508 
152000. 1259, 
167000. 1272. 
44 BREAKPOINTS SUPERHEATER 
0.0 0.0 
20000. О 
28000, 0 
34000. 20 
390000 4.0 
45000, 5.2 
48000, 6.0 
50000. 6.5 
52000. 7.0 
54000, f 
56000, 8,0 
8000, 8.6 
60000. 9.2 
62000. 10.0 
68000. 10.5 
660006 11.1 
68000. 12, 0 
70000 12.8 
74000 14.0 
76000. 150 
78000, 15.9 
80000. 16.6 
84000. 18.0 
86000» 19.0 
88000. 20.0 
90000. 21.0 
94000» 23.0 
96000. 28.0 
98000. 2540 
100000. 26.0 
108000. 30.0 
117000, 25.0 
120000. 37.0 
125000. ‚ 4040 
130000. 44.0 
140000, 50.0 
145000. 54.0 
150000. 59.0 
160000» 66.0 
165000» 70.0 
170000. 74.0 
180000. 84.0 
190000. 92.0 
200000. 101.0 
=01 
50 152000. 
DOl TIME INCREMENT 
55.0 TOTAL TIME 
le 
0.2 Ut PLOT TIME 
0194 VARIBLES TO BE PRINTED 


O30BRIEN GH DELTA SYSTEM 
AIR FLOW(C 94) (PERCENT) 


9994 


VS TIME 
147 


(SEC) 


CURVE 2 


(90 РЕВ) (МЕС КАМР ) 4-4 





(LBS/HR) 


L κ 28) 
С-89) ιδ) ον ἽΙ 


LTA FEEDWATER FLOW 
LTA FUEL FLOW (C-8 


шш 


D 
D 





Appendix IV 


Subroutines INVAL and DIAGRAM and Data for Cruising 
Conditions for Ten to Ninety Per Cent Full Power 
Simulation of the Complete Boiler 





t 


SUBROUTINE INVAL 
DIMENSION G(99)5H(99)5.»P(99),5Q( 99) 9R(99).C( 99) ,CO( 99) 5 COO( 99) 
1С000(949)»С0000(99)»Т(10)»Х1(900)»Х2(900)»Х3(900)» 
2Ү1(900) »Ү2(900)»ҮЗ(900)» 
SITITL& 12) sJTITLETT2) KT ITLE O2 
4Х(1540)»Ү(15»40) 
СОММОМ |6»Н»Р»0»Б»С>С0»С00»С000»С0000»7»Х1»Х2»Х3»Х»Ү»Ү1»Ү2»ҮЗ 
І 


С(50) TS READ FROM IC CARD 
C(42)215900, 

C(97)=150 

С(92)=С(50) 


CALL CURVE (2.47.50) 

CALL CURVE (1576550) 

CALL ADDER (4859489763239-47) 
С(77)-0.04%С(48)-22.0 
C(49)=C177) 
C(52)=C197)-C(49) 
C(53)2(G(53)/H( 53) ) *C( 52) 
С(78)-0.0001142 *C(50)+3.0 
С(57)=С(78) 
C(54)=C(53)+C(57) 
Е(55)-С(54) 

C(67)2C(55) 
C(66)2C(67)-3.0 
C(3)=C(66)*(H(66)/G( 66) ) 
C(65)=C(3) 

С(79)-С(65) 

ШІ 75)-С(79)/,254 
ЕГТІ2ІЕС(75)/6(75) 
C(56)=C(55)=-C(67) 

CALL CURVE(4, 73, 72) 

CALL CURVE (9588572) 

CALL CURVE (5574,88) 


C(71)=C(72) 

C(70)=C(71)/C(73) 

C(69)=C(70) ; 
C(68)=C(69) 

CALL CURVE(6»63»68) 

C(61)=C(63) 

C(60)2C(61) 


E 581 -CI561*(GL581/P(581] 
C(51)2C( 58) *G(l51] 

C(59)» Ct60)-C( 51] 
IF(C(55)-Ct67))15,152 
C(64)=C(55) 

GO TO 3 

C(64)=C(67) 

CONTINUE 

FUEL OIL FLOW (GF) FOR CRUISING SPEED =4140. 
C(^6)2C(64)-34.0 > 
C(45)=C (46) #H(46) /G(46) 
С(39)<С(45)/0.00141 
C(41)=C(42)-G(40)*C(39) 
C(40)2C(39)*G( 40) 


C(90)2C(39Y 
C(38)02C(^41)-C( 39] 
C(87)2C(38) 


C(81)30400141*C( 41) *340 150 





C(80)= С(38)%0.00141 +3.0 
C(44)=C(81) 


C(43)=C(80) 
C(5)=C(1%4)+15 e 0 
C(33)=C(5) 

C(93)=C(50) 
C(85)=C(93)+C(28) 
С(4)-0.000114%С(93)43.0 
C(29)=C(&) 

C(30)=C (29) 

43 Wie C57 ) =C (30) 
C(32)=C(31)%*(G(32)/H(32)) 
C(2&)=C(97)+C(32)—C(33) 
С(25)-С(24)%(6(25)/Р(25)) 
КЕТУКМ 

ЕМО 


15 





1000 
2000 


3000 
4000 


1100 
2100 


3100 
4100 


1200 
2200 


3200 
4200 


SUBROUTINE DIAGRAM 


ОІМЕМБТОМ 6(99)»Нн(99)»Р(99)»0(99)»ЕК(99)»>С(99)»С0(99)»С00(99)» 
1С000(99)»С0000(99) ЭТ(10)»Х1(900) »Х2(900) ,Х3(900)» 
Ἐν ο ο αι 

BITITLE(12) ώΤΙΤΙΕ 11) Ακ τει 
4X%(15940) s¥(15940) 

COMMON GsHoP 9QsRsCsC0sCO00sCO00 » C0000 »T oX19X29X39X9Vo9VYlsV29Y3 

IF (C(99)-23,0) 10510520 

C(50) 516600, 4(5915.0%С(99)) bee 

CONTINUE | 

CALL CURVE (2947550) 

CALL CURVE (1976550) 

CALL ADDER (48,48 5769239-47) 

CO(77)=C(77) 
C(77)=e04*C (48) - 33.0 
CALL CMPLXPL(49949977) 
С(97)-15. 
CALL ADDER (523525979~4950) 
CALL REALPL (535953952) 
СО(78)-С(78) 

С(78)-.0001142%С(50) 43.0 

CALL CMPLXPL(57»57»78) 

CALL ADDER (5455495355750) 

CALL CMPLXPL (55955954) 

CALL ADDER (565565555-6750) 

CALL CMPLXPL (58,58,56) 

CALL CURVE (988,72) 

CALL GAIN (51951558) 

IF(C(51)-30,0 )1000»1000»2000 

IF (C(51)-04e0) 300053000 34000 

С( 51) =30,0 

GO TO 4000 

С(51)-0.0 

CONTINUE 

CALL REALPL (59959958) 

IF (C(59)—-3020) 1100911002100 

IF (C(59)-020) 310093100 94100 

С(59)-20.0 3 

GO TO 4100 

С(59)-0.0 

CONTINUE 

CALL ADDER(6O» 605 5195990) 

IF (C(60)-30¢): 1200120092200 
ІЕ(С(60)-0.) 3200932004200 
C(60)=308 

GO TO 4200 

C(60)=08 

CONTINUE 

CALL CMPLXPL (61»61»60) 

CALL DERIVA (62362961) 

CALL ADDER (63563961 s6230) 

THE AIR SYSTEM GOES HERE 

CALL CURVE (3568963) 
CALL REALPL (69569568) 
CALL CMPLXPL (70570969) 


G(71)=C (73) 

CALL GAIN: (71971970) 
G(72)=C(74)" 
H(72)=C(74) 


“CALL REALPL( 72972971) 


СА. СОВУМЕ (4›73›72) 
CALL CURVE (5,74,88) 
CALL GAIN (75975972) 

С0( 79.) =С (79) 

С (79) =С (75) *.254 
С(94&)=С(75)*#1,2 
CALL REALPL(65»65,79) 

CALL DELAY1 (353,65) 
CALL REALPL (6656653) 
C(67)2Cl66) 43,0 
ІҒ(С(55)-С(67))191»2 

1 С(64)-С(55) 
GO TO 3 

2 C(64)=C(67) 

3 CONTINUE ` 

C OIL SYSTEM GOES HERE 

КОО) 
C(34)=C(64)—-(C( 46) 43-0) 
CALL REALPL (35535534) 
CALL GAIN (36936934) 
CALL ADDER (379379353 3630} 

CALL REALPL (86986937) 
CALL ADDER(389385—86 98730) 
CALL ADDER( 39939-38941 50) 
CALL ADDER(893893399~9090) 
CALL GAIN(40»40939) 

CALL ADDER(41»415-4054250) 

СО(81)-С(81) 

С0(80)-С(80) | 
С(803- С(38)%0.00141 42.0 
С(81)-0.00141%С(41)%2.0 
CALL CMPLXPL(43543580) 
CALL CMPLXPL(44 944981) 
CALL ADDER (455945944 5-—-4390) 
CALL REALPL (46946545) 

C WATER SYSTEM 

(0 5)=C (5) 
C(5)=Ct14)+15.0 
CALL CMPLXPL (3393395) 
CALL ADDER (245245-33932997) 
CALL CMPLXPL (84984924) 
CALL GAIN(25925 584) 

CALL REALPL (26526: 84) 
CALL ADDER (2792732532690) 
CALL DELAY1 (252527) 
CALL CMPLXPL(2852852) 
C(85)2C(93)4C( 28) 
socia) 
C(4)=02000114*(C(93)4+C(28) ) +340 
CALL CMPLXPL (2992994) 
CALL REALPL (30530929) 
CALL ADDER(315»315575—-3050) 
CALL REALPL (32932531) 

С WATER LEVEL 
CALL ADDER(91»91550,-9250) 
CALL REALPL (696991) 
CALL REALPL( 797991) 
CALL CMPLXPL-( 89897} 
CALL ADDER (9999897630) 
CALL REALPL (10510528) 





С 


CALL CMPLXPL(11»11»10) 
CALL REALPL(12512589) 
CALL СМРЕХРЕ (13.13.12) 
CALL DELAY1 (1515123) 

CALL ADOER (145145951151) 


STEAM PRESS 


CALL CMPLXPL (15.15.01 | 
CALL REALPL (16516915) 

CALL ADDER (1791751591690) 
CALL REALPL (825982589) 

CALL REALPL (185918582) 
CALL REALPL (194919918) 

CALL ADDER (209209189~19390) 
CALL REALPL (21521528) 

CALL REALPL (22.22.21) 
CALL ADDER (235,23,205»-175»-22) 
С(95)-С(39)/13.1 
C(96)=(Ct93)+Ct28))/1660% 
RETURN 

END 

END 


s4 





С DATA FOR 10-90 PER CENT FULL POWER IN 23 SECONDS 


-48 NUMBER OF BLOCK DATA CARDS 


О] „ 36 
02 0.25 
03 0.045 


06 1202 -06 


07 2.805 -05 0.3141 
08 1.0 0.763 
10 5.54- 07 
id eo 0.422 
12 1.86 -04 2.22 
ο O 0.76 
15 2.000368 1.4545 - 
16 „028 
18 7.67 -03 .0445 
19 1.732 1.76 
21 1„73 -06 
22 1.0 0.0445 
25 1.26 
26 0.1 
28 9360. .531 
29 19.2 2.62 
20. 0,222 0.222 
32 147 29.4 
33 4.0 2.4 
25 0.292 
36 2.82 
40 0.0445 
43 30.3 5.9 
44 30.3 5.9 
46 30.87 29.4 
49 50,0 „б 
5] .944 
D3 73.6 2904 
25 250. 95 
ЕТ” 1942 2.62 
58050333, 13% 2 
59 106 
61 100. 7.9 
62 0289 
65 2.5 2.5 
66 30.87 29.4 
69 1.0 1.0 
το 2.0 0.56 
72 000865 ο 15775 
75 1.134 
82 1.0 
83 1.0 
84 1.25 0.675 
86 24804. 31.2 
09 
10 ВКЕАК РТ5 
0.0 
16600. 
34000. 
56000. ғ 
66500. 
100000. 


1.588 
0.495 


O o 4 


©0627 


. 781 
1902 


250% 


19.2 
333,03 


100. 


INI 5 


BOILER PRESSURE 
1200. 


1201. 


1203. 
1208. 
1211. 
1226. 
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CURVE 


1 





120000. 1237, 
180000. 1250, 
1520006 1259, 
167000, 2725 
40 PRESSURE DROP ACROSS THE SUPERHEATER 
0.0 0.0 
16600. | ὃ 
28000. 2.0 
34000, 3,0 
39000, 4.0 
45000» 5.2 
48000, 6.0 
50000. 6.5 
52000. 7.0 
54000 7.4 
56000, 8,0 
60000, 902 
62000, 10.0 
640006 10.5 
66000. 11.1 
68000, 12.0 
70000, 12, 8 
74000 14.0 
76000» 15.0 
78000.» 15.9 
80000.» 16.6 
84000. 18.0 
900006 21.0 
940006 23.0 
98000. 25.0 
100000» 26.0 
108000» 30.0 
117000. 35.0 
120000. 37.0 
125000. 40.0 
130000, 48.0 
140000, 50.0 
145000. 54.0 
152000» 59.0 
160000» 66.0 
165000. 70.0 
170000, 74.0 
180000, 8440 
190000. 92.0 
200000. 101.0 
10 BREAKPONTS FOR CURVE 3 
-100. 1440. 
240 1440. 
5.0 1440,0 
6.0 1600.0 
10.25 2800.0 
14.3 4300,60 
18.5 6400.40 
22.75 9200.0 
27.0 12800.0 
1000» 12800, 
08 BLOWER 0 
6.2 00436 
5522 000436 
| ‚00762 


20.8 
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CURVE 4 





9:2 «00735 
63.2 . «00538 
78.8 00488 
100. е 00488 
Олык SPEED VS TIME CONSTANT (1000 RPM) CURNE 5 
9 οὐ 4 
1.0 ο0 4 
2.0 \ 00585 
3.0 4 . 09 
4.0 „128 
5.0 E 9162 
6.0 4214 
. ο „272 
8.0 „3015 
20.0 3015 
10 
0.0 3.0 
1440. 3.0 
1440. 5.0 
1600. 6.0 
2800. 10.25 
%300. 14.3 
6400. 18.5 
9200. 2 2. ub 
12800, 27.0 
18000. 31.3 
08 GPR VS М 
0.0 2.73 
3.0 2813 
be 2.23 
De 1.77 
6. 1.26 
Ге 1,0 
δ. 1.0 
10. 1.0 
08 INTEGRAL GAIN VS N 
0.0 06498 , 
3.0 2498 
4.0 6498 
5.0 2498 
6.0 .397 
7.0 2397 
8.0 2444 
10.0 . 444 
08 BLOWER SPEED VS ОВ CURVE 10 
6.2 1.8 
5.35 1.8 
20.8 2.912 
34.7 4.14 
49.2 5.107 
63.2 6.708 
74.8 7.480 
100. 7.480 
-01 
50 16600. 
ο ο. τ 
50.0 TOTAL TIME 
1. 
oll 
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0196 

03 J FENICK BOILER 10-90 WITH CR 
9994 PERCENT OF AIR 

9995 PERCENT OF FUEL 

9996 PER CENT OF WATER 
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